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ABSTRACT 


This  thesis  involves  investigation  of  linear  filtering 
models  as  a  means  of  generating  texture  in  images.  Various 
autoregressive  filter  models  are  used  to  generate  various 
textures,  and  the  results  are  analyzed  to  determine  rela¬ 
tionships  between  filter  parameters  and  texture  characteris¬ 
tics.  A  two-dimensional  counterpart  to  the  autoregressive 
integrated  moving  average  (ARIMA)  model  from  one-dimensional 
time  series  analysis  theory  is  developed  and  tested  for 
texture  modeling  applications.  All  these  models  are  driven 
by  white  noise,  and  to  the  extent  that  real  images  can  be 
reproduced  this  way,  advantages  in  image  texture  transmis¬ 
sion  could  be  realized.  Results  of  this  work  indicate  that 
the  purely  autoregressive  models  work  well  for  some  types  of 
image  textures,  but  that  for  the  textures  studied  the  ARIMA 
model  is  not  particularly  suitable. 
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I.  INTRODUCTION 


The  purpose  of  this  thesis  is  to  investigate  the  types 
and  quantity  of  image  textures  generated  using  a  two- 
dimensional  (2-D)  extension  of  the  Autoregressive  Integrated 
Moving  Average  (ARIMA)  model.  For  the  one-dimensional 
(e.g.,  time  series)  case,  the  theories  and  formulas  describ¬ 
ing  this  model  are  outlined  in  Box  and  Jenkins  [Ref.  l:pp. 
85-103].  The  1-D  ARIMA  model  is  useful  when  the  time  series 
to  be  modeled  is  not  stationary  but  exhibits  some 
homogeneity  in  the  sense  that,  except  for  statistical 
differences  between  parts  of  the  time  series,  these 
different  parts  of  the  process  behave  similarly.  In  these 
cases  some  suitable  difference  of  the  process  may  be 
stationary,  and  hence  may  be  accurately  modeled  by  an 
Autoregressive  Moving  Average  (ARMA)  or  a  purely 
Autoregressive  (AR)  model.  The  resulting  stationary  time 
series  (generated  by  an  appropriate  ARMA  or  AR  filter  with 
white  noise  input)  is  applied  to  an  integra-tion  or 
summation  filter  (the  inverse  of  the  difference  operation) 
to  generate  the  original  nonstationary  time  series.  [Ref. 
l:p.  85]  Figure  1-1  shows  a  block  diagram  of  this  process. 

This  work  attempts  to  extend  these  concepts  to  two- 
dimensional  signal  processing.  In  order  to  simplify  the 
model,  the  moving  average  (MA)  portion  of  it  will  be 
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Figure  1-1  Block  Diagram  for  the  Autoregressive 
Integrated  Moving  Average  Model 


eliminated  (i.e.,  no  zeros  in  the  filter  Z  transform)  so 
that  only  purely  AR  models  will  be  considered  for  stationary 
image  generation.  The  procedures  for  modeling  image 
textures  using  AR  models  with  white  noise  input  are  well 
established  [Ref.  2:pp.  454-456].  However,  a  suitable  two- 
dimensional  difference  operation  and  its  inverse  must  be 
found  to  implement  the  concepts  outlined  above. 

The  research  is  divided  into  four  areas: 

1)  Investigation  of  the  various  types  of  image  textures 
generated  using  AR  models  where  filter  coefficients 
and  size  are  determined  a)  arbitrarily,  b)  using  a 
two-pole  separable  model,  and  c)  using  a  four-pole 
separable  model.  Separability  refers  to  the  fact  that 
the  Z  transform  of  the  AR  filter  can  be  factored  into 
components  representing  each  dimension  or  direction  of 
the  image. 

2)  Selection  of  a  difference  operator  and  a  realizable 
inverse  (integration  or  summation)  filter. 

3)  Application  of  the  above  autoregressively  generated 
images  to  the  summation  filter,  and  evaluation  of 
these  results. 

4)  Attempted  reproduction  of  actual  images  textures  using 
AR  models  whose  coefficients  are  determined  using  the 
statistics  of  the  image,  and  comparison  of  these 
results  to  those  obtained  by  a)  applying  the 
difference  operator  to  the  real  image,  b)  finding  the 
coefficients  of  the  AR  model  that  reproduce  the 
difference  image,  and  c)  applying  the  difference  image 
to  the  summation  filter.  This  comparison  was  intended 
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to  discover  what  improvements,  if  any,  may  be  realized 
using  the  ARIMA  model  vice  a  purely  ARMA  (or  AR) 
model . 

The  remainder  of  the  thesis  is  organized  as  follows. 
Chapter  II  contains  methods  and  results  of  investigating 
various  autoregressive  image  models.  Chapter  III  deals  with 
the  formulation,  development,  and  testing  of  the  two- 
dimensional  summation  filter.  Chapter  IV  contains  the 
results  of  applying  various  AR-generated  images  to  the 
summation  filter  and  addresses  the  application  of  the  ARIMA 
model  to  real  image  textures.  Chapter  V  outlines  conclu¬ 
sions  on  the  results  and  applicability  of  the  ARIMA  model. 
Although  the  ARIMA  modeling  was  not  highly  successful  in 
reproducing  the  textures  studied  here,  plausible  reasons  are 
given  for  their  failure  and  conjectures  are  made  about  those 
circumstances  where  the  model  would  be  more  successful. 
Appendix  A  provides  information  on  the  computer  algorithms 
used  to  implement  the  equations  governing  the  above 
processes.  Appendices  B  through  G  contain  derivations  of 
spectral  and  autocorrelation  equations,  and  the  correspond¬ 
ing  graphical  results,  governing  the  AR  processes  in  Chapter 
II.  Appendices  H  and  I  contain  graphical  results  associated 
with  the  inverse  filter  development  in  Chapter  III. 

Image  data  were  generated  using  computer  programs 
written  in  FORTRAN,  compiled  using  Version  4.5  under  the 
VAX/VMS  Version  4.4  operating  system.  The  images  were 
displayed  on  the  COMTAL  Vision  0ne/20.  The  gray  level 
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II.  THE  AUTOREGRESSIVE  IMAGE  MODEL 


A  two-dimensional  signal  (such  as  an  image  texture)  can 
be  modeled  using  a  two-dimensional  AR  model  with  white  noise 
input.  The  governing  equations  in  the  spatial  domain  are  of 
the  following  form  [Ref.  3:pp.  325-326]: 

P-1  Q-l 

y(n,m)  =  -  £  a-j^y  (n-i,m-j )  +  w(nfm)  (2.1) 

i=0  j=0 
(i ,  j  ) ^  (0 , 0) 

where  y(n,m)  is  a  signal  representing  the  generated  image 
texture  at  pixel  location  (n,m) ,  a^j  is  the  filter  coeffi¬ 
cient  matrix,  and  w(n,m)  is  a  two-dimensional  white  noise 
signal.  The  system  function  corresponding  to  the  filter  of 
Eq.  (2.1)  is  given  by 


Y(zlfz2)  = 


-1 


-1 


-1-1 . 


1+a10Zl  +a01Z2  +allZl  Z2  +"'+aP-lQ-lZl 


-(P-1)  -(Q-l) 


W(zlfz2) 


(2.2) 


where  and  z2  are  the  Z  transform  variables  corresponding 
to  spatial  coordinates  n  and  m.  Ideal  white  noise  has  an 
autocorrelation  function  that  is  an  impulse  and  a  flat  (con¬ 
stant)  power  spectrum  with  magnitude  corresponding  to  the 
variance  of  the  white  noise  process  [Ref.  4:pp.  22-26]. 
Therefore  determination  of  the  filter  coefficients  a^j  will 
define  the  generated  image  process.  Procedures  will  be 
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outlined  later  to  estimate  the  coefficients  from  real  image 
data.  At  this  point  analytical  methods  will  be  used  to 
select  these  coefficients  and  the  resulting  images  will  be 
studied.  Figure  2-1  shows  an  example  of  a  white  noise  input 
image . 


A.  ARBITRARILY  SELECTED  FILTER  COEFFICIENTS;  P  =  2,  Q  =  2 
In  order  to  get  an  initial  idea  of  what  types  of  images 
might  be  generated  using  a  2  x  2  AR  filter  with  white  noise 
input,  filter  coefficients  were  at  first  selected  arbitrari¬ 
ly,  but  subject  to  a  stability  constraint.  The  primary  con¬ 
straint  on  coefficient  selection  is  that  of  filter 
stability.  Using  the  DeCarlo-Strintzis  Theorem  dealing  with 
multidimensional  filter  stability  [Ref.  3:pp.  197-198], 

alternately  setting  =  1  and  z2  =  1  and  determining  the 
location  of  the  pole  in  the  remaining  dimension  will 
indicate  whether  or  not  the  filter  is  stable.  If  the 
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magnitude  of  the  pole  in  the  remaining  dimension  is  less 
than  1,  the  filter  is  stable.  Even  with  this  condition, 
however,  there  are  an  infinite  number  of  possible  filter 
coefficient  combinations.  The  additional  constraint  of 
a10  =  aoi  can  be  used,  and  comparisons  of  results  using 
various  values  of  a1;L  can  be  made. 

Figure  2-2  shows  the  form  and  directionality  convention 
used  for  the  autoregressive  filter,  along  with  the  corres¬ 
ponding  difference  equation  and  its  Z  transform. 


n(z2) 


h(n,m) 


f 

a01 

all 

1 

a10 

CT  — 

nfZj,) 


1  1 

y  (n,m)  =  -  [  l  a. . -y(n-i,m-j)  +  w(n,m) 

i=0  j=0  13 
(i,j)^(0,0) 


Y(z^,Z2)  H(z^,Z2) •W(z1 ,z0)  — 


\'“2'  ~  -1  -1  -1  -1 
1+a10Zl  +a01Z2  +allZl  Z2 


W(z1,z2) 


Figure  2-2  Autoregressive  Filter  Impulse  Response, 
Difference  Equation,  and  Z  Transform 


Although  it  is  difficult  to  make  precise  predictions  in 
two  dimensions,  one  can  expect  that  the  sign  and  magnitude 
of  a^o  or  a0i  would  influence  the  correlation  between  pixels 
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in  the  corresponding  directions.  For  example,  a  positive 
value  for  a10  ®ight  be  expected  to  yield  an  image  with 
substantial  variation  in  the  n  direction  (low  correlation) , 
particularly  if  the  magnitude  of  a^o  is  near  1.  A  negative 
value  for  aQi  with  magnitude  near  1  might  yield  an  image 
with  lower  variations  in  pixel  intensity  (high  correlation) 
in  the  m  direction.  Since  the  filter  is  not  necessarily 
separable  (i.e.,  the  denominator  of  H(z1,z2)  cannot  be 
factored  into  the  form  D(z^) • D(z2) ) ,  conclusions  drawn  from 
this  line  of  reasoning  may  not  be  completely  correct. 

Initial  attempts  at  generating  images  with  arbitrarily 
selected  coefficients  yielded  rather  uninteresting  results 
having  very  little  contrast  or  discernible  pattern. 
Continued  experimentation  with  combinations  where 
a10  m  a01  <  0  eventually  yielded  more  interesting  image  tex¬ 
tures.  Figures  2-3  and  2-4  show  the  results  of  using  the 
constraint  a^j  =  a01  =  -0.35  and  various  values  of  a^  for 
the  filter  coefficients.  For  positive  values  of  a13_,  the 
images  are  rather  "grainy,”  with  higher  magnitudes  yielding 
a  somewhat  "finer"  graininess.  There  are  also  some  overall 
intensity  differences  observed.  For  the  negative  values  of 
aH,  the  results  are  much  more  interesting.  As  the  magni¬ 
tude  increases,  there  is  a  gradually  more  noticeable  upper 
left  to  lower  right  orientation  of  the  image  texture,  and 
the  variations  from  lower  left  to  upper  right  become 
smoother  as  well.  Using  a^g  =aoi  =  -0.38  and  a^  =  -0.24 
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yields  minor  variations  in  texture  pattern  and  overall  image 
intensity  when  compared  to  the  previous  case. 

Using  a^Q  -  a01  »  0.35  and  varying  a^  from  0.0  to  0.8 
(Figure  2-5),  the  images  obtained  deviate  very  little  from 
the  mean  intensity  value,  and  possess  minor  differences  in 
graininess.  With  these  positive  coefficients  some  negative 
correlation  might  be  expected,  and  the  fact  that  these 
images  are  "grainy"  indicates  the  existence  of  some  negative 
correlation  or  high  spatial  frequency  characteristics.  On 
initial  examination,  however,  the  low  contrast  of  the 
generated  images  tends  to  obscure  the  observed  graininess. 

B.  TWO  POLE,  SEPARABLE  AUTOREGRESSIVE  MODEL 

In  general,  it  is  difficult  to  relate  the  nature  or 
properties  of  a  two-dimensional  filter  to  the  precise  nature 
of  an  image  texture  that  may  be  generated  when  white  noise 
is  applied  to  that  filter.  In  order  to  simplify  the  effort 
and  to  obtain  a  better  understanding  of  the  problem,  the 
case  where  the  filter  (and  resulting  image  texture)  are 
separable  is  considered.  For  the  two  pole  separable  case 
considered  here,  the  filter  transfer  function  can  be 
factored  into  expressions  in  z^  alone  and  z2  alone.  The 
expressions  in  z^  and  z2  each  have  one  pole  on  the  real  axis 
in  their  respective  Z  domains.  Figure  2-6  illustrates  the 
filter  structure,  the  corresponding  difference  equation,  and 
its  Z  transform. 
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Figure  2-5  Images  Generated  Using  Arbitrarily 
Selected  Filter  Coefficients 
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Va'a'a'a'a ->'a  aV-Vj. 


n(*l) 


i=0  j-0 
(i,j)*(0,0) 


<*11  *  *10'*01> 


Y(z1,z2)  =  H(zirz2) •W(z1#z2)  = 


1«1021  1+a01z2 


•  W(z  ,z2) 
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Figure  2-6  Autoregressive  Filter  Form,  Difference 
Equation,  and  Z  Transform 

Here  it  is  relatively  easy  to  relate  stability  of  the 
filter  to  the  location  of  the  poles  in  the  z^  and  z2  planes. 
Since  the  quarter  plane  filter  is  separable  and  the 
components  are  causal,  one-dimensional  filter  stability 
theory  can  be  used  to  state  that  the  poles  in  each  plane 
must  have  magnitude  less  than  1  to  ensure  filter  stability. 
Figures  2-7  through  2-9  show  images  resulting  from  this 
model  for  various  values  of  a^o  and  a01-  Note  that  the  sign 


Figure  2-9  Images  Generated  Using  a  Two-Pole 
Autoregressive  Model 
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Careful  comparison  of  the  images  resulting  from  various  com¬ 
binations  of  a^o  and  a01  leads  to  the  following 
observations: 

1)  When  poles  are  located  in  the  same  place  in  the  z-±  and 
z2  planes  on  the  negative  side  of  the  real  axis, 
magnitudes  near  1  yield  a  fine  graininess  with  patchy 
areas  and  low  overall  contrast.  As  the  magnitude  of 
the  pole  decreases,  the  graininess  becomes  more  coarse 
and  the  result  is  more  like  the  original  white  noise 
input.  No  directional  quality  in  the  image  pattern  is 
observed. 

2)  When  poles  are  located  in  the  same  place  on  the  posi¬ 
tive  side  of  the  z^  and  z2  real  axes,  a  somewhat  dif¬ 
ferent  result  is  observed.  For  magnitudes  near  1,  an 
image  of  patchy  light  and  dark  areas  results,  with 
differing  amounts  of  correlation  between  pixels  in 
different  areas.  Slightly  discernible  "lines”  in  both 
the  horizontal  and  vertical  directions  are  also 
observed.  For  lower  pole  magnitudes  on  the  positive 
side  of  the  real  axis,  the  decreased  effect  of  the 
filter  on  the  white  noise  input  is  again  observed. 
This  result  is  more  like  the  white  noise  and  has  more 
contrast  than  the  corresponding  result  using  poles  on 
the  negative  side  of  the  real  axis. 

3)  For  pole  placements  in  the  z1  and  z2  planes  which  are 
on  the  negative  side  of  the  real  axis  and  are  of 
unequal  magnitude,  the  results  have  a  very  fine 
graininess  and  low  contrast.  Some  slight  direction¬ 
ality  is  observable  in  the  image  patterns,  with  lower 
frequency  variations  evident  in  the  direction  corres¬ 
ponding  to  the  pole  with  smaller  magnitude. 

4)  For  pole  placements  in  the  z^  and  z2  planes  which  are 
on  the  positive  side  of  the  real  axis  and  are  of 
unequal  magnitude,  much  more  directionality  and  varia¬ 
tion  is  observable  in  the  image  pattern. 

5)  As  the  poles  are  placed  on  opposite  sides  of  the  real 
axis  and  are  separated  by  a  greater  distance,  direc¬ 
tionality  becomes  more  evident  (with  higher  frequency 
variations  in  the  direction  of  the  more  negative 
pole) .  As  the  pole  separation  becomes  greater  and  as 
the  pole  magnitudes  become  closer  to  l,  smoother  sinu¬ 
soidal  variations  are  evident. 

6)  When  the  pole  magnitudes  are  equal  and  have  opposite 
sign,  the  image  generated  using  pole  magnitudes  close 
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to  1  exhibits  high  frequency  sinusoidal  variations  in 
the  direction  of  the  negative  pole.  The  image  gener¬ 
ated  with  the  lower  magnitude  poles,  as  would  be 
expected  from  the  above  results,  resembled  the 
unfiltered  white  noise. 

In  order  to  explain  these  image  patterns  analytically, 
analysis  of  the  power  spectrum  and  autocorrelation  function 
of  this  process  is  useful.  Since  this  model  is  separable, 
the  analysis  can  be  conducted  in  each  direction  separately. 
The  power  spectrum  is  defined  by  [Ref.  4:pp.  24-34]: 

Sy(w)  -  a2|H(e>))|2  =  a  2H(e  ^H(e-3W)  (2.3a) 
where 


H(e>)  =  H ( z )  |  z=e?u  (2.3b) 

and  for  this  case 

H(z)  =  - i—r  (2.3c) 

1  +CIZ 

Here  a2  is  the  magnitude  of  the  white  noise  power  spectrum. 
We  can  assume  that  a2  =  1  with  no  loss  of  generality  of  the 
results,  since  o2  does  not  affect  the  shape  of  the  frequency 
response . 

The  autocorrelation  function  is  related  to  the  filter 
transfer  function  through  the  equations  [Ref.  5:pp.  391- 

395) : 
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h(n) 

=  Z_1[H(z) ] 

(2.4a) 

y(n) 

=  h(n)*w(n) 

00 

(2.4b) 

>i 

=  a  2  l  h (n) *h (n-£ ) 

(2.4c) 

n=-°° 


Specific  forms  of  the  power  spectrum  and  autocorrelation 
function  are  given  in  Appendix  B.  Since  Ry(£  )  =  Ry(~M 
[Ref.  5 :p.  388],  calculating  the  expression  for  Ry(£),  £  <  0 

is  not  necessary.  From  Appendix  B,  the  results  are: 

Sy  (w )  =  J  (2.5) 

1  +2aoos(aj)  +a 

Ry(£)  =  -^2-  I  >  o  (2.6) 

1  -a 

Appendix  C  shows  the  results  of  these  equations 
graphically  for  various  values  of  a  .  The  relationship  be¬ 
tween  the  power  spectra  and  their  corresponding  autocorrela¬ 
tion  functions  conforms  to  the  expected  results  from  theory 
(i.e.,  low  frequency  spectrum  with  smooth  autocorrelation 
function,  and  high  frequency  power  spectrum  with  rapidly 
varying  autocorrelation  function)  [Ref.  6:pp.  139-142].  The 
plots  in  Appendix  C  also  demonstrate  that:  1)  For  poles  on 
the  positive  side  of  the  real  axis  in  the  Z  plane  low  fre¬ 
quencies  predominate  and  for  poles  on  the  negative  side  of 
the  real  axis  high  frequencies  predominate,  and  2)  Lower 
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magnitudes  of  a  result  in  a  broader  power  spectrum  and  a 
wider  range  of  frequencies  of  significant  magnitude.  Both 
of  these  observations  agree  with  the  image  results.  For 
images  generated  using  a  more  negative  pole  in  a  given 
direction,  fine,  high  frequency  graininess  is  observed  in 
that  direction  (though  the  low  contrast  or  low  variation 
about  the  mean  intensity  may  tend  to  make  this  effect  less 
noticeable) .  When  a  more  positive  pole  is  used,  lower  fre¬ 
quency  variations  are  more  evident  in  the  corresponding 
direction.  As  lower  magnitudes  are  used  for  a  in  a  given 
direction,  more  random  variations  (indicative  of  a  wider 
range  of  significant  frequency  components)  are  observed  in 
that  direction.  The  form  of  the  autocorrelation  function 
for  these  cases  approaches  the  autocorrelation  function  for 
white  noise,  i.e.,  an  impulse.  Negative  poles  should  yield 
high  frequencies  since  the  negative  side  of  the  real  axis  in 
the  Z  plane  represents  a  digital  spatial  frequency  of  tt  , 
while  positive  poles  in  the  Z  plane  correspond  to  a  digital 
spatial  frequency  of  zero.  Note  that  even  when  the  poles 
are  placed  such  that  a  spatial  frequency  of  zero  should  pre¬ 
dominate,  there  are  some  low  frequency  random  variations  in 
the  resulting  images.  Since  the  power  spectra  of  the  posi¬ 
tive  poles  all  contain  some  non-zero  frequency  components 
(they  are  not  perfect  impulses  at  zero) ,  this  characteristic 
is  expected. 
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C.  FOUR  POLE,  SEPARABLE  AUTOREGRESSIVE  MODEL 


For  the  cases  considered  in  this  section,  H(z)  again  can 
be  factored  into  expressions  in  z^  and  Z2-  However  here 
each  factor  is  a  2nd  degree  polynomial  with  two  poles  in  the 
denominator.  Figure  2-10  illustrates  the  filter  structure, 
the  applicable  difference  equation,  and  the  corresponding  Z 
transform. 


n(z2) 
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Figure  2-10  Autoregressive  Filter  Form,  Difference 
Equation,  and  Z  Transform 


Since  all  of  the  a^j  coefficients  are  real,  the  poles 
must  1)  both  be  on  the  real  axis,  or  2)  occur  in  complex 
conjugate  pairs  in  the  and  z2  planes.  Again,  pole  magni¬ 
tudes  must  be  less  than  1  to  ensure  filter  stability.  We 
will  assume  here  that  the  poles  in  each  of  the  factors  have 
equal  magnitudes  and  opposite  (or  0  or  ± tt)  phase.  Letting 
a i  =  magnitude  of  poles  in  the  z^  plane 

0 i  =  pole  angle  (phase)  in  the  z^  plane 

a  2  =  magnitude  of  poles  in  the  z2  plane 

0  2  =  pole  angle  (phase)  in  the  z2  plane 

and  using  Euler's  relation,  the  denominators  of  H(Z!)  and 
H(z2)  can  be  expressed  as  follows: 


-1  -2 

l+allZl  +a20Zl  =  (1-aie 


j®l  _i  ~-^l  -1  -1  2  -2 

z1  )  (1-a.je  z1  )  =  l-2a^cos (6^)  mz^  -Kx^ 


1+a01Z2  +a02Z2  (1_a2e 


j®2  -1  -1  -1  2  -2 

z^)  (l-ct2e  z^)  ~  1_2a2ac>s(0o)z'>  +a-z' 


j2'2  22 


Hence : 

a10  =  ”20t1*cos(8  j)  a01  =  -2a2 •  cos (  02) 

2  2 
a20  “  al  a02  "  a2 

This  gives  a  relationship  between  the  pole  magnitude  and 


angle  in  the  Z  domain  and  the  filter  coefficients  in  the 
spatial  domain. 


1.  Complex  Conjugate  Poles 

Figures  2-11  through  2-13  illustrate  images 
generated  using  this  model  for  various  complex  conjugate 
pole  combinations  in  the  z^  and  z2  directions.  In  comparing 
each  of  these  image  textures  in  terms  of  the  relative  effect 
of  pole  positioning  in  each  direction,  the  following  obser¬ 
vations  can  be  made: 

1)  For  images  generated  using  poles  of  equal  magnitude 
and  angle  in  both  directions,  graininess  with  no 
directionality  to  the  pattern  resulted.  Higher  pole 
angles  yielded  finer  (higher  frequency)  graininess  and 
less  contrast.  Lower  magnitude  poles  yielded  a  more 
random  and  less  structured  graininess  pattern  at  the 
same  pole  angle. 

2)  Using  a  pole  angle  of  zero  (pole  on  positive  real 
axis)  in  one  direction  and  a  pole  of  some  non-zero 
angle  in  the  other  direction  yielded  images  with 
highly  directional  sinusoidal  patterns.  The  direction 
corresponding  to  the  pole  on  the  real  axis  was  not 
totally  devoid  of  variation,  but  variations  were  slow, 
i.e.,  of  very  low  frequency.  The  spatial  frequency  of 
the  sinusoidal  pattern  can  be  increased  by  increasing 
the  pole  angle.  Large  magnitude,  high  pole  angle 
combinations  yielded  much  cleaner  and  more  structured 
textures  than  low  magnitude,  low  pole  angle  combina¬ 
tions.  Lower  magnitude,  high  pole  angle  combinations 
yielded  less  structured  textures  where  directionality 
was  evident  but  the  sinusoidal  pattern  was  obscured. 
Low  magnitude,  low  pole  angle  combinations  yielded 
very  random,  non-structured  textures  of  relatively 
high  contrast. 

3)  Using  poles  in  the  z^  and  z2  planes  with  the  same 

magnitude  but  different  pole  angles  resulted  in  some 
directionality  if  there  was  a  sufficiently  large 
magnitude  and  difference  in  the  pole  angles.  As 

observed  earlier,  the  direction  with  the  higher  pole 
angle  had  the  higher  spatial  frequency.  Large  pole 
magnitudes  (close  to  1)  resulted  in  more  structured 
but  rather  low-contrast  images  (the  low  contrast 
seemed  to  obscure  the  high  frequency  nature  of  the 
pattern  somewhat) .  Pole  angles  in  z^  and  z2  that  were 
of  close  value  made  it  difficult  to  detect  the  higher 
frequency  (higher  pole  angle)  direction.  Low 
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magnitudes  basically  negated  the  pole  angle  effects 
and  yielded  a  very  random,  unstructured,  high  contrast 
texture. 

4)  Where  pole  magnitudes  were  close  and  pole  angles  were 
different  in  the  z^  and  z2  planes,  some  directionality 
in  the  texture  was  observed.  Again,  the  high  pole 
angle  direction  yielded  the  highest  frequency.  When 
the  pole  angles  in  both  Z  domains  had  similar  values 
and  the  magnitudes  of  the  poles  differed,  graininess 
with  little  or  no  discernible  directionality  or 
structure  resulted. 

These  observations  are  consistent  with  classical 
pole-zero  frequency  response  analysis  [Ref.  7:pp.  323-331]. 
There  is  a  direct  relationship  between  pole  angle  and 
spatial  frequency  in  a  given  direction,  and  the  magnitude  of 
the  poles  affects  the  amount  of  structure  and  definition  of 
the  sinusoidal  pattern  of  a  given  frequency  in  a  given 
direction.  Higher  magnitude  poles  result  in  a  narrower 
bandwidth  of  the  filter  and  yield  more  structure  and 
sinusoidal  pattern  definition.  Low  pole  magnitudes  give  the 
filter  wider  bandwidth  and  yield  images  with  less  structure 
and  definition  and  more  randomness  in  a  given  direction. 
While  directional  dependencies  are  evident  given  pole 
magnitude  and  angle  in  a  given  direction,  it  does  not  appear 
that  a  pattern  in  one  direction  is  totally  independent  of  a 
pattern  in  the  other  direction.  This  would  be  expected, 
even  though  the  model  is  separable,  due  to  the  cross  terms 
in  the  filter  structure. 

Filter  power  spectrum  and  autocorrelation  analysis 
can  be  conducted  in  this  case,  as  in  the  case  of  the  two 
pole  model.  The  derivations  for  Sy(^)  and  Ry(.)  are 
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somewhat  more  involved,  and  are  given  in  Appendix  D.  The 
resulting  expressions  for  Sy(^)  and  Ry(  i)  from  Appendix  D 
are: 


Sy(^) 


(2.7) 


Ry(M 


a 2  ,cosU9)  cos(  (2+2)9)-a  cos(£9)  x 

- y-( - 5 - 3 - T  ’ 

2sin  9  1  -a  1+a  -2a  cos(26) 


(?-  >  0) 


(2.8) 


The  plots  of  these  functions  for  the  various  pole 
magnitudes  and  angles  used  are  given  in  Appendix  E.  The 
9=0  case  is  equivalent  to  having  2  poles  on  the  real  axis 
at  a  given  magnitude  in  the  Z  plane.  As  would  be  expected, 
the  power  spectrum  for  each  model  showed  higher  magnitudes 
at  digital  frequencies  close  to  the  pole  angle.  Higher  pole 
magnitudes  yielded  sharper,  more  well-defined  power  spectrum 
magnitudes  at  the  given  frequency,  and  lower  pole  magnitudes 
yielded  less  well  defined  more  spread-out  power  spectra. 
Low  pole  magnitudes  almost  completely  obliterated  evidence 
of  low  frequency  power  spectrum  components,  and  degraded  its 
definition  and  sharpness  at  higher  frequencies.  This 
corresponds  to  the  observed  results  in  the  image  textures 
generated.  The  autocorrelation  functions  also  reflect  the 
appropriate  relationship  to  the  power  spectra  as  outlined  in 
the  discussion  for  the  two  pole  case,  i.e.,  greater 
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variation  in  the  autocorrelation  function  indicates  greater 
variation  between  pixels  a  given  distance  apart,  which  in 
turn  implies  higher  spatial  frequencies. 

2.  Two  Real  Poles 

Rather  than  using  complex  conjugate  pole  locations 
to  obtain  real  filter  coefficient  values,  two  poles  on  the 
real  axis  may  also  be  used  for  a  given  direction.  They  may 
be  placed  at  different  locations  on  the  real  axis,  or  they 
may  be  placed  together.  The  latter  situation  is  equivalent 
to  the  0=0  (or  6  =  ±n  if  placed  on  the  negative  real  axis) 
case,  as  mentioned  above, 
relevant  equations  are: 


(l-aaz_1) (l-e^z-1) 


where 


For  the  two  real  pole  case,  the 


: - =i - 

1+a01Z  +a02Z 


a01  =  "  (aa  +  lb) 

a02  =  1  a 'a  b 

For  these  experiments  a  transfer  function  of  the 
complex  conjugate  pole  form  was  used  for  the  direction, 

with  =  0.9  and  !  =  0.  For  the  z2  direction  a  transfer 
function  with  two  poles  on  the  real  axis  was  used.  The 


image  textures  that  result  for  various  values  of  a  a  and  ^ 
are  given  in  Figures  2-14  and  2-15.  Some  observations  can 
be  made  about  these  results: 

1)  With  poles  placed  at  the  same  value  on  the  Z2  real 
axis,  rather  unstructured,  low  frequency  variations 
are  observed  in  the  Z2  direction.  The  more  positive 
poles  result  in  very  slow  variation  in  the  image 
texture,  while  the  lower  magnitude  positive  poles  show 
more  variation  in  the  z2  direction.  When  the  poles 
are  placed  in  the  same  location  on  the  negative  side 
of  the  real  axis,  a  low  contrast  image  with  some 
noticeable  high  frequency  variations  results. 

2)  As  the  poles  are  moved  farther  apart  on  the  Z2  real 
axis,  high  frequency  variations  with  increasing  struc¬ 
ture  and  oscillatory  form  are  evidenced  in  the  z2 
direction. 

3)  When  poles  with  equal  magnitude  and  opposite  sign  are 
used,  fairly  structured  high  frequency  variations  are 
evidenced  in  certain  areas  of  the  image,  while  low 
frequency  variations  are  evident  in  other  areas  in  the 
z^  direction.  Higher  magnitude  poles  yield  more 
discernible,  structured  variations,  while  lower 
magnitude  poles  of  opposite  sign  yield  discernible  but 
non- oscillatory  high  frequency  variations  in  certain 
areas  of  the  image. 

Of  particular  interest  is  the  fact  that  two  poles 
placed  at  the  same  value  on  the  negative  real  axis  in  the  Z2 
plane  yielded  some  high  frequency  variations.  This  is  in 
keeping  with  the  fact  that  values  on  the  negative  real  axis 
correspond  to  a  pole  angle  (and  corresponding  digital 
frequency)  of  8  =  tt  .  The  presence  of  poles  on  the  negative 
side  of  the  real  axis  of  the  Z2  F  ne  seems  to  give  rise  to 
the  high  frequency  variations  with  gradually  more  structure 
and  oscillatory  appearance  as  the  pole  is  moved  to  the  left 
(more  negative) . 
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Figure  2-15  Images  Generated  Using  a  Four  Pole 

(Two  Real  Poles)  Autoregressive  Model 
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The  expressions  for  the  power  spectrum  and  autocor¬ 
relation  function  for  the  random  process  produced  by  driving 
the  filter  of  Figure  2-10  with  white  noise  are  derived  in 
Appendix  F.  It  is  shown  there  that  the  power  spectral 
density  and  the  autocorrelation  function  are  given  by: 


Sy  (tO  ) 


Ry(M 


2  2  2 
1-2  ^aa"l^+'aaat>+aaab ) 003  +2aaa^cos  (2oo)  +a^ 

+  2VVWo 


t+l  ,  t+1 

aa  ^3% 

^“b 


2+t 

A* 


(2.9) 

(2.10) 


Plots  of  these  functions  for  the  various  values  used 
in  this  section  are  given  in  Appendix  G.  The  power  spectrum 
results  are  consistent  with  the  observed  image  spatial  fre¬ 
quency  characteristics.  Both  low  and  high  frequency 
components  were  contained  in  some  of  the  power  spectra,  and 
were  manifested  in  the  corresponding  images  as  both  low  and 
high  frequency  variations  in  the  z2  direction.  The  nature 
of  the  autocorrelation  functions  related  to  the  power 
spectra  that  contained  low  and  high  frequency  components  was 
interesting.  Autocorrelation  functions  with  much  variation 
but  all  positive  values,  rather  than  the  equal  magnitude 


positive  and  negative  values  evidenced  in  earlier  results, 
seems  to  reflect  the  higher  level  of  correlation  related  to 
the  low  frequency  (smoother  variations)  aspect  of  the  image 
texture  variations. 


D.  IMAGE  TEXTURE  ROTATION  TRANSFORMATION 

If  an  image  signal  x(n1,n2)  consists  of  a  rotated 
version  of  another  image  w(m2,m2)  such  that  n^  =  Im^  +  Jm2 
and  n2  =  Km^  +  Lm2,  where  I,  J,  K,  and  L  are  integers  and 
IL-KJ  7*  0,  then  the  Z  transform  X(zlfz2)  is  given  by 
W(z£, z2, z^, z2)  [Ref.  3:p.  182].  A  45°  rotation  corresponds 
to  1  =  1,  K  =  1,  J  =  1,  L  s  -i,  If  we  use  the  four  pole 
separable  result  for  H(Zi,z2),  as  shown  in  Figure  2-10,  and 
apply  the  above  transformation  (z^  z^z^,  z2  ■+•  z^z21)  ,  we 
find  after  simplification: 


Hp  ( Z  iZ  o )  —  zii  HI  Ho  2  — l  — I  ^2  —3  — 1 

1+a01zl  z2  +a02Zl  Z2+a10Zl  Z2  +a10a01Zl  +a10a02zl  Z2 


“2  2  “"3  — 1  _.r  ^ 

+  a20zl  Z2  +a^.0a01Zl  Z2  +a20a02zl 


(2.11) 


Notice  that  this  transfer  function  is  not  separable  but 
consists  of  a  rotated  version  of  a  separable  filter.  Figure 
2-16  illustrates  the  support  of  the  denominator  polynomial 
for  this  filter.  It  has  the  form  of  a  non-symmetric  half¬ 
plane  infinite  impulse  response  (HR)  filter,  so  it  is 
recursively  computable. 
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Figure  2-16  Rotation  Transformation  Filter  Form 


The  application  of  this  filter,  using  filter  coefficients  of 

■h71 

-D*> 

the  four  pole  separable  filter  with  poles  at  -*■  0.9e  , 

%2  O.Se^0  in  the  original  separable  filter  yielded  the 
result  shown  in  Figure  2-17. 


Figure  2-17  Result  of  Rotation  Transformation 


E.  SUMMARY 

Autoregressive  models  can  produce  a  variety  of  image 
textures.  For  general  two-dimensional  models,  the  system 
functions  are  generally  not  factorable  and  singularities 
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occur  on  surfaces,  not  at  isolated  points.  For  these 
reasons  it  is  difficult  to  design  two-dimensional  filters 
for  images  and  predict  the  resulting  character  of  the 
images.  Indeed,  even  to  ensure  stability  of  the  filter  is 
not  trivial.  As  a  result  we  concentrated  here  on  separable 
forms,  which  by  their  nature  are  much  easier  to  analyze. 
Certain  types  of  texture  patterns  using  various  separable 
autoregressive  models  can  be  predicted  based  on  filter  pole 
placement  in  the  and  z2  planes.  Arbitrary  or  random 
selection  of  filter  coefficients  can  yield  interesting  but 
generally  unpredictable  results.  Obviously,  an  infinite 
number  of  variations  on  the  models  above  could  be  attempted. 
Ultimately,  the  anticipated  utility  of  the  textures 
generated  will  guide  the  process  of  model  and  parameter 
selection. 


IN  TWO  DIMENSIONS 

To  implement  the  2-D  AJRIMA  model,  the  inverse  of  a 
filter  representing  a  suitable  difference  operator  is 
needed.  One  possible  2-D  difference  operator  is  the 
Laplacian,  which  has  the  impulse  response  shown  in  Figure  3- 


1  [Ref.  8 : pp .  212-213]. 

m 

h(n,m) 

0  -1 

0 

— 

-i  0  -i 

0  -1 

0 

Figure  3-1  Laplacian  Impulse  Response 

Its  implementation  involves  convolving  it  with  an  image 
and  is  represented  by  the  following  difference  equation: 

1  1 

y(n,m)  =  l  l  bi^x(n-i,  m-j )  (3.1) 

i=-l  j=-l 


where 


x(n,m)  is  the  image  input  signal,  and 

b-H  is  the  filter  coefficient  matrix 
(bij  =  h (i , j ) ) 
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In  the  Z-transform  domain  this  can  be  written  as: 


Y(z1#z2)  =  H(z1( z2) -X(zlf z2) 

=  (4  -  z"1  -  zx  -  z~l  -  z2)*X(z1,z2) 


In  areas  of  an  image  where  adjacent  pixels  have  similar  gray 
levels  (low  frequency,  homogeneous  areas),  the  result  of 
this  operator  will  be  approximately  zero.  Where  significant 
or  sharp  differences  in  gray  levels  between  adjacent  pixels 
exist,  the  result  of  this  operation  will  be  farther  from 
zero.  Thus  the  Laplacian  difference  operator  is  sometimes 
used  as  an  "edge  detector." 

The  problem  addressed  in  this  chapter  is  constructing 
the  inverse  of  the  operator.  In  the  Z  domain,  the 
expression  for  the  inverse  would  be  [Ref.  4:p.  36]: 


H-1(zi,z2) 


1 

n(zltz2) 


which  has  an  expansion  as  an  infinite  series  of  positive  and 
negative  powers  of  z^  and  z2.  That  is,  considering  this 
expression  as  a  problem  in  long  division,  the  result  of  such 
division  would  be  an  expression  of  the  form: 


H"1(z1,z2) 


CO  CO 


where 


a^j  *  coefficient  values  of  z^fi  resulting  from  the 
long  division  t 


Note  that  if 


1 


4  -  z'1  -  Zi  -  z-1  -  z2 


)  a.  z  iz  J 

1  ]  1  2 

1  =  -*-  3  =  -,: 


then  cross  multiplication  would  yield: 


(4 


—  7~ i  -  7 i  —  7~ ^  — 


Zl 


Z2)  \  ' 


1  =  -'”  ]  =  -• » 


-1  -1 
af  j  z  z 


1  (3.3) 


The  double  summation  expression  in  z^  and  z2  will  be 
truncated  and  considered  to  be  an  FIR  filter  with  finite 
support  and  coefficients  a^j .  This  approximates  the  desired 
inverse  filter.  In  particular,  we  will  use  the  following 
constraints: 

1)  Choose  the  limits  of  summation  to  be  equal  in  both 
directions,  i.e., 


-j 

2 


This  results  in  a  "square"  region  of  support  for  the 
filter  (all  values  outside  assumed  zero) . 

2)  Force  the  values  for  the  filter  coefficients  to  be 
symmetr ic,  i.e.,  a  ^  ^  a_  ^  j  3 i— j  a  _  ^  ^  a  j  ^ 

=  a_ji  =  a j =  a_ j . 

Using  these  constraints  and  implementing  the  cross 
multiplication  equation  (3.3)  will  result  in  an  expression 
in  z^  and  z2,  with  each  combination  of  the  z^z"^  terms  having 
a  coefficient  whose  form  is  a  summation  of  terms  in  a^j 
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The 


where  the  coefficients  of  a^j  are  either  4  or  -1. 
coefficient  of  the  z^z®  term  must  equal  1  and  the  coeffi¬ 
cient  of  any  other  z^xz~p  term  must  equal  zero  to  satisfy 
equation  (3.3)  . 

As  an  example,  let  L  =  1.  Equation  (3.3)  can  then  be 
expressed  as: 

11 

(4  -  z-1  -  Zi  -  z'1  -  z 2 )  ’  l  1  a  j.  j  z -]1  z ^ 

i=-l  j=-l  J 

1  1  11 
=  l  +  crzi  +  0‘z2  +  o-z-j^  +  ... 


Performing  the  double  summation  yields 


i=-l  j=-1 


a^.^z^zj-  +  a_ioziz2  +  a-nzfz--1  + 


'1  2 


+  a00z°z0  +  ao-i.z°z^  +  ai.iz-lzl  +  a^zf1 


+  a01Z2 1  +  allzllz21 


Performing  the  cross  multiplication  would  yield  45  different 
terms  in  various  combinations  of  z^z^.  Combining  these 
terms  to  find  the  coefficient  expression  for  each  z^xz^p 
term  soon  becomes  rather  tedious  and  impractical  for  even 
moderate  values  of  L.  An  alternative  way  to  proceed  is  to 
choose  a  zaz^  term  on  the  right  hand  side  of  the  equation, 
and  for  each  term  in  the  expression  4-z'j1-z1-z  21-z2  >  deter¬ 
mine  what  values  of  i  and  j  are  required  so  that  when 
each  term  is  multiplied  by  a^jz'^-z”},  it  will  result  in  an 
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expression  in  the  chosen  zaz^  term  on  the  right  hand  side  of 
the  equation.  Choosing  z°z°  ^  on  the  ri<3ht  hand  side  of 
the  equation,  we  have: 


4 '  ai  j 

z-iz-j 
z  i  z2 

=  CZ^Z1^ 

when 

i 

=  0 

and 

j 

=  0 ; 

so 

o 

o 

it! 

Ti¬ 

ll 

u 

z~xz~i 

1  2 

=  cz°z° 

1  2 

when 

i 

=  - 

1  and 

j 

=  0; 

so 

c  =  -a_ 

-Z1  *  ai  j 

z ""  i  z  “  D 
Z1  Z2 

=  cz°z° 

1  2 

when 

i 

=  1 

and 

j 

=  0; 

so 

c  =  -a 10 

z"1  •  a  i  t 

2 

=  cz°z° 

1  2 

when 

i 

=  0 

and 

j 

=  -l  ; 

so 

c  =  _a0-l 

~z2 ‘ ai j 

Zliz23 

=  cz°z° 

1  2 

when 

i 

=  0 

and 

j 

=  1 ; 

so 

c  =  -a01 

So  the  coefficient  for  z°z°  is  simply  a  summation  of  the  c 

12 

terms  obtained  above,  i.e., 


(4a00~  3.10  ~  a10  “  a0-l  “  a0l)z°z2 

This  entire  expression  must  equal  1  to  satisfy  (3.3),  and 
since  z°z°  =  1,  4a0o  ~  a-io  "  a10  “  a0-l  ~  a01  =  1  also- 

Using  the  same  method  for  the  z^z^  term  on  the  right 
hand  side  yields: 


4*aijz-iz-3 

=  czj-z^ 

when 

i 

=  -1  and  j  =  -1 ;  so  c  = 

4a_i_i 

z-1’  ai-iZ-:Lz_3 

1  J  1  2 

=  CZ^-Z1 

1  2 

when 

i 

=  -2  and  j  =  -1;  so  c 

=  -a_2 

-Zi  '  &\AZ~XZ~') 

1  1  2 

=  cz^z1 

1  2 

when 

i 

=  0  and  j  =  -1 ;  so  c  = 

-a0_i 

Z2  1’aijzilz  ~2 

=  cz^z1 

1  2 

when 

i 

=  -1  and  j  =  -2;  so  c 

=  _a-l 

~z2’aijz~1izp 

=  «l2»2 

when 

i 

=  -1  and  j  =0;  so  c  = 

:  _a-10 
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The  resultina  term  in  z}z\  is: 

1  z 

( 4a-l-l  *  a-2  - 1  ~  a0-l  "  a  - 1  -  2  '  a-10)zl1ZT 

This  expression  must  equal  zero,  since  there  is  no  z1z1  term 
on  the  right  side  of  (3.3),  so  4a_1_1  -  a_2-i  -  a0-i  -  a-l-2 
~  a-10  =  °* 

This  procedure  can  be  extended  to  any  number  of  z~,1z~J 
terms.  When  this  is  done,  the  resulting  expressions  for  the 
coefficients  of  z^lz"3  can  be  formed  into  a  set  of  simultan¬ 
eous  equations  in  order  to  solve  for  the  a^j  coefficient 
values.  However,  due  to  the  symmetry  condition  imposed 
above,  some  of  the  equations  for  the  coefficients  of  the 

terms  are  linearly  dependent.  For  values  of  l  and  j 
that  yield  unique  or  distinct  values  for  a^j,  the  resulting 
z'^z'j  coefficient  expressions  are  linearly  independent. 
For  example,  the  coefficient  expression  for  the  z'^z®  term 
is  linearly  independent  of  the  coefficient  expression  for 
the  z~xz~®  term,  since  a10  *  an-  But  the  coefficient  ex- 
pression  for  the  z^z°2  term  is  linearly  dependent  on  the 
coefficient  expression  for  the  z^z"1  terTn>  since  a10  =  aoi- 
Using  only  the  linearly  independent  equations  for  a  given 
filter  size  yields  a  set  of  p  equations  in  p  unknowns,  where 
p  is  the  number  of  unique  and  distinct  a^j  values  in  the  in¬ 
verse  filter.  The  value  of  p  is  related  to  the  size  of  the 
desired  inverse  filter.  If  the  size  of  the  filter  is  N  •  N, 


the  number  of  unique  a^j  values  using  the  symmetry  of 
constraint  above  is: 


=  (~  +1)  +  (~)  +  -1) 


...  +  1 


(3.4a) 


(N+l ) (N+3) 

4 


(3.4b) 


For  example,  with  N  =  7,  the  unique  a^j  values  in  a  7-7 
inverse  filter  can  be  represented  by  a33,  a 32 ,  a31,  a3(3, 

a22 i  a21'  a2 0 '  a 1 1 >  a i o ,  a00.  Though  there  are  49  elements 
in  a  7-7  filter,  all  of  them  are  equal  to  one  of  these 
values  listed,  due  to  symmetry.  Obviously,  N  is  constrained 
to  be  odd  for  a  square  filter  with  a  unique  element  a0Q  in 
the  center. 

Tne  solution  of  the  resulting  p  equations  yields  the 
values  for  the  p  filter  elements  or  coefficients.  This 
defines  the  FIR  approximation  to  the  inverse  filter.  It  is 
only  an  approximation  due  to  the  finite  size  constraint 
imposed,  and  it  might  be  expected  that  the  larger  the  filter 
size,  the  better  the  approximation. 

An  algorithmic  procedure  for  obtaining  the  a^j 
coefficients  is  outlined  below.  An  example  follows. 

1)  Determine  the  desired  size  of  the  inverse  filter. 

2)  For  each  combination  of  (positive)  i,j  values  corres¬ 
ponding  to  a  unique  a^j  filter  coefficient,  identify 
the  five  term  summation  equation  associated  with  each 
z~llzp  term. 
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3)  Combine  equal  a^j  values  and  develop  a  matrix  of 
coefficients  for  the  a^j  values.  Let  this  matrix  be 
A. 

4)  Denoting  the  vector  of  unique  a^j  values  as  a,  the  set 
of  simultaneous  equations  in  matrix  form  is: 


where 

a  =  (a00  a10  axl  a2o  *21  a22  •••  ^-l  n-1 

~2~  ~T 

The  top  row  of  A  corresponds  to  the  summation  of  a^  j 
terms  that  represents  the  coefficient  of  the  zvz^ 
term. 

5)  Solve  (3.5)  for  a. 

An  example  of  this  procedure  is  appropriate  at  this 
point.  For  an  inverse  filter  of  size  Nx  N: 

Step  1 

Let  N  =  5  (therefore  L  =  2) 

Step  2 

The  coefficients  corresponding  to  each  unique  z-1z"J 
term  are: 

z°z°  -*  4a00  -  a10  -  a_10  -  a0i  -  *o-l 

z-ij-O  ^  4a10  -  a20  -  a00  -  axl  -  a1.1 

z^z-1  -*■  4alx  -  a2i  -  a01  -  a12  “  a10 
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Z~2Z° 

zi  Z2 

4a20 

~  a30 

~  a10 

"  a2 1  ”  a2-l 

r  2Z-1 

1  2 

-  4a21 

"  a3 1 

“  all 

"  a2 0  "  a2-2 

,-2--2 

'1  Z  2 

-*■  4a22 

”  a32 

“  a12 

“  a23  -  a2 1 

Step  3 

Combining  equal  terms  in  Step  2  and  expressing  the 
coefficients  of  a^j  in  matrix  form  yields  an  A  matrix  of: 


Step  4 

With  A  given  in  Step  3,  the  a  vector  for  (3.5)  is 

a  =  Ca00  a10  all  a2 0  a2 1  a22  J 

Step  5 

Solving  (3.5)  for  a  involves  inverting  A  and  multiplying 
it  by  [1  0  0  0  0  0]T.  Thus: 

a  =  A-1  *  [  1  0  0  0  0  0]T 
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Appendix  H  illustrates  the  forms  of  the  resulting 
inverse  filters  of  various  sizes,  as  well  as  the  normalized 
and  unnormalized  filter  cross  sections. 

One  way  to  validate  the  resulting  inverse  filter  is  to 
convolve  it  with  the  original  Laplacian  difference  operator. 
The  result  should  approximate  an  impulse  at  the  origin. 
Appendix  I  shows  the  results  of  this  convolution  using  3^3, 
5x5,  7  »  7,  9  x  9,  15  x  15,  and  21  x  21  size  inverse  filters. 
It  is  seen  there  that  as  the  size  of  the  filter  gets  larger, 
it  becomes  a  better  approximation  to  the  true  inverse  and 
the  convolution  looks  more  like  an  impulse. 

To  test  the  application  of  this  filter  on  an  actual 
image,  a  test  image  was  filtered  using  the  Laplacian 
difference  operator.  Then  the  resulting  image  data  were 
filtered  again  using  various  size  inverse  filters.  The 
results  are  shown  in  Figures  3-2  and  3-3.  Note  that  the 
image  resulting  from  Laplacian  FIR  filtering  seems  more 
stationary  than  the  test  image,  which  was  one  of  the  desired 
results.  Inverse  filtering  of  that  result  yields  images 
that  are  progressively  more  similar  to  the  original  test 
image  as  the  size  of  the  inverse  filter  increases.  However, 
a  rather  large  inverse  filter  is  needed  to  accurately 
reproduce  the  image.  The  result  of  the  21x21  inverse 
filter  is  quite  similar  to  the  test  image,  with  some  loss  of 
contrast  or  darkness  in  certain  areas,  but  with  essentially 
the  same  pattern.  The  effect  of  the  size  limitation  of  the 
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7*7  inverse  filter  9x9  inverse  filter 


15  *  15  inverse  filter 


21  x  21  inverse  filter 


Figure  3-3  Results  of  Filtering  Test  Image  with 
Laplacian  FIR  Filter  and  its  Inverse 
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inverse  filter,  as  manifested  in  the  convolution  of  the 
Laplacian  and  its  inverse  above,  would  seem  to  explain  the 
lack  of  perfect  test  image  reproduction.  Larger  inverse 
filter  sizes  could  be  tried,  but  large  inverse  filter  sizes 
relative  to  image  size  would  result  in  a  significant  portion 
around  the  edge  of  the  image  having  only  a  part  of  the 
filter  applied  to  it.  This  would  adversely  affect  the 
overall  quality  of  image  reproduction. 
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IV.  APPLICATION  OF  THE  ARIMA  MODEL  TO  IMAGE  TEXTURES 

As  outlined  in  Chapter  I,  the  utility  of  the  ARIMA  model 
centers  around  the  fact  that  a  difference  operator  applied 
to  an  image  texture  may  improve  the  stationarity  of  the 
image  statistical  characteristics.  A  stationary  image 
texture  is  required  for  accurate  modeling  by  autoregressive 
techniques,  and  it  was  hoped  that  application  of  the 
autoregressively  generated  texture  to  an  approximate  inverse 
of  the  difference  operator  may  yield  a  more  accurate  or 
recognizable  representation  of  the  original  nonstationary 
image,  as  compared  to  a  purely  autoregressively  generated 
version. 

A.  APPLICATION  OF  LAPLACIAN  INVERSE  FILTER  TO  AUTO¬ 
REGRESSIVELY  GENERATED  IMAGES 

As  an  initial  examination  of  the  effects  of  the  inverse 
filter  developed  in  Chapter  III  on  image  textures,  selected 
images  generated  in  Chapter  II  were  input  to  the  21*21 
version  of  that  filter.  Figures  4-1  through  4-3  illustrate 
the  results.  All  attempts  resulted  in  a  blurred  or  smoothed 
version  of  the  original  image.  Since  the  inverse  of  a 
difference  operation  is  a  summation  or  "integration" 
operation,  and  since  integration  operations  can  be  expected 
to  smooth  or  blur  (low  pass  filter)  signals  [Ref.  8:pp.  136- 
154],  the  results  are  not  surprising.  However, 
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except  to  the  extent  that  blurring  is  useful  or  desirable, 
applying  the  summation  filter  to  image  signals  that  are  not 
based  on  the  application  of  the  corresponding  difference 
filter  to  that  signal  seems  to  be  of  little  utility. 

In  the  remainder  of  this  chapter  we  consider  applica¬ 
tion  of  the  summation  filter  to  regenerate  actual  image 
textures . 


B.  AUTOREGRESSIVE  FILTER  PARAMETER  ESTIMATION  PROCEDURES 

The  first  step  in  testing  the  ARIMA  model  is  to  estimate 
the  autoregressive,  quarter  plane  filter  parameters  required 
to  model  the  real  image  textures  and  the  signal  resulting 
from  application  of  the  Laplacian  operator  to  those  images. 
For  a  zero-mean  signal,  these  model  parameters  are  found  by 
solving  a  set  of  Normal  equations.  In  these  equations  the 
white  noise  covariance  is  referred  to  as  the  prediction 
error  covariance.  The  Normal  equations  can  be  expressed  as 


R 


(  4.1) 


ap-l 


0 


where  the  R  matrix  is  the  correlation  matrix  for  the  signal 
(block  Toeplitz  with  Toeplitz  blocks) ,  the  a  vector  consists 
of  appropriately  ordered  filter  coefficient  vectors,  and  s 
is  a  vector  containing  the  prediction  error  covariance  as 


the  first  and  only  nonzero  element.  Here  =  [a^g  aii 
ai2  .  .  .  a^  q-i]T  and  T  =  [a  2  0  0  ...  0]T. 

Calculating  the  correlation  matrix  and  prediction  error 
covariance  from  the  image  signals,  and  solving  (4.1)  for  the 
~5T  vectors,  provides  all  the  parameters  needed  for  the  2-D  AR 
model.  The  multichannel  form  of  the  Levinson  recursion  can 
be  used  to  solve  these  equations  more  efficiently  [Ref.  2:p. 
454]  . 

C.  APPLICATION  TO  REAL  IMAGE  TEXTURES 

Actual  image  textures  used  here  are  from  the  image  data 
base  at  the  University  of  Southern  California's  Signal  and 
Image  Processing  Institute  [Ref.  9:pp.  13-14].  The  images 
selected  from  this  data  base  are  contained  in  a  book  by 
Brodatz  [Ref.  10].  Portions  of  the  images  of  size  128  x  128 
pixels  were  obtained  and  used  as  a  basis  for  processing. 
Filter  coefficients  were  calculated  for  the  real  image 
textures  shown  in  Figures  4-4  and  4-5.  Various  filter  sizes 
were  tried  to  determine  which  yielded  the  best  results  in 
generating  a  particular  image,  and  a  quarter-plane  filter 
size  of  4  x4  was  selected.  Results  of  autoregressive 
filtering  of  white  noise  using  the  appropriate  calculated 
coefficients  to  model  each  texture  are  given  in  Figures  4-6 
and  4-7.  Images  generated  by  applying  the  Laplacian  differ¬ 
ence  operator  to  the  real  images  are  shown  in  Figures  4-8 
and  4-9.  Autoregressive  generation  of  these  images  using 


Gravel 


Grass 


Wall 


After  Laplacian  Filtering 


filters  with  the  corresponding  calculated  coefficients  are 
given  in  Figures  4-10  and  4-11.  Finally,  the  application  of 
the  signal  represented  by  the  images  in  Figures  4-10  and  4- 
11  (without  the  0-255  scaling  reflected  in  these  figures)  to 
the  21  x  21  inverse  filter  described  in  Chapter  III  yields 
the  images  shown  in  Figures  4-12  and  4-13.  Comparison  of 
all  of  the  above  results  yields  the  following  observations: 

1)  Autoregressive  modeling  of  the  water,  grass  and  sand 
textures  yielded  good  results.  Some  of  the  other  tex¬ 
tures  with  more  structure  and  sharp  local  variations 
were  not  reproduced  well. 

2)  Autoregressive  reproduction  of  images  created  after 
application  of  the  difference  operator,  with  the 
exception  of  the  water  image,  yielded  generally  poor 
results.  As  observed  in  Chapter  III,  the  application 
of  the  difference  operator  produces  a  seemingly  more 
stationary  result,  but  the  edge  structure  that  remain¬ 
ed  in  most  of  the  images  after  application  of  the  dif¬ 
ference  operator  was  in  general  not  reproducible  using 
a  purely  AR  model. 

3)  Application  of  the  inverse  filter  to  the  image  signal 
generated  by  AR  model  reproduction  of  the  difference 
operator  results  yielded  smoothed  versions  of  those 
results.  This  is  similar  to  what  was  observed  in  Sec¬ 
tion  A  of  this  chapter  when  images  were  applied  to  the 
inverse  filter  that  were  not  based  on  the  specific 
data  generated  by  the  difference  operator. 

As  a  final  test  of  the  ARIMA  model,  a  64  x  64  contrast 
enhanced  aerial  photograph  of  trees,  with  smoother 
variations  and  in  general  less  edge  structure  than  the  other 
images  tested,  was  tried.  The  results  are  shown  in  Figure 


4-14.  Though  this  image  seemed  somewhat  better  adapted  to 


the  model,  overall  the  same  observations  outlined  above 


apply. 
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Figure  4-11  Laplacian  Filtered  Image  Textures 
Generated  Using  an  AR  Model 
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Figure  4-12  Image  Textures  Generated  Using  an 
ARIMA  Model 


Grass 


Figure  4-13  Image  Textures  Generated  Using  an 
ARIMA  Model 
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Figure  4-14 


Final  Test  of  ARIMA  Model  on  Contrast 
Enhanced  Trees  (Magnification  X2) 


D.  SUMMARY 


The  effectiveness  of  AR  reproduction  of  image  data  using 
white  noise  input  and  filter  coefficients  calculated  based 
on  the  statistics  of  the  image  signal  is  highly  dependent  on 
the  nature  of  the  image  data.  The  water  image,  for  example, 
with  its  smoothly  varying  and  rather  homogeneous  nature,  was 
quite  well  adapted  to  AR  reproduction.  Other  images  with 
more  structure,  abrupt  variations,  and  more  non-homogeneous 
characteristics,  were  not  autoregressively  reproducible  to 
any  great  extent. 

Using  the  ARIMA  model,  it  seems  that  the  operation  of 
the  inverse  filter  is  very  sensitive  to  the  nature  of  the 
input  data.  Input  data  that  are  strictly  based  on  the 
difference  operator  output  can  reproduce  the  original  image, 
as  was  found  in  Chapter  III.  However,  the  AR  model  used  to 
generate  the  inverse  filter  input  (based  on  the  statistics 
of  image  signal  produced  using  the  difference  operator)  does 
not  generate  image  data  accurately  enough  to  reproduce 
images  that  resemble  the  real  images  tested. 


V.  CONCLUSIONS 


This  thesis  sought  to  explore  experimentally  and  to 
understand  how  linear  filtering  models  could  be  used  to 
generate  texture  in  images.  Of  particular  interest  was  the 
investigation  of  2-D  ARIMA  models  to  see  if  they  might  be  of 
any  utility  in  this  effort.  Some  time  was  spent  exploring 
separable  2-D  models  to  understand  how  transfer  function 
pole  placement  affected  image  texture  characteristics. 
Image  textures  generated  using  these  models  and  applied  to 
the  summation  filter  yielded  blurred  or  smooth  textures  with 
seemingly  little  variety  or  utility.  The  ultimate  test  of 
the  model  was  its  ability  to  reproduce  actual  image 
textures.  The  purely  AR  portion  of  the  model  reproduced  a 
few  types  of  actual  textures  well.  However,  the  full  ARIMA 
model  failed  to  generate  image  textures  that  resembled  the 
source  images  used.  Many  of  the  textures  had  strong  edge 
differences  that  were  not  accurately  reproducible  by  the  AR 
model.  Also,  the  summation  filter  developed  seemed  very 
sensitive  to  deviations  in  image  signal  data  from  that 
generated  by  application  of  the  difference  operator;  that 
is,  the  procedure  seemed  not  to  be  "robust." 

Since  many  of  the  images  tested  here  have  definite  edge 


structure , 

the  difference 

image  had 

lines  which  were 

not 

reproduced 

well  by  the 

AR  mode  1  . 

Correspondingly 

the 

integrated  AR  model  did  not  reproduce  the  original  image. 
For  this  type  of  image,  a  combination  of  a  line  point 
process  model  [Ref.  11]  with  the  integrator,  would  possibly 
have  been  more  suitable.  The  image  of  trees  had  not  such 
edge  structure  and  produced  somewhat  better  results. 
Further  experimentation  with  images  of  this  type  and  the 
ARIMA  model  would  perhaps  be  appropriate. 


APPENDIX  A 


COMPUTER  PROGRAMS.  SUBROUTINES.  AND  FUNCTIONS 

Listed  below  are  the  names,  associated  computer  systems, 
and  functions  of  the  various  computer  algorithms  used  to 
accomplish  this  thesis  research.  All  programs  were  written 
by  the  thesis  author  unless  otherwise  noted.  Program  source 
codes  are  given  at  the  end  of  this  appendix  (except  for  the 
MAKFI L*  series) . 

A.  PROGRAMS 

1.  AL'TOREG  (VAX/ VMS  FORTRAN) 

The  program  did  the  following: 

1)  Generated  a  128  ■  128  zero  mean  white  noise  matrix- 
using  subroutine  PGAUSS. 

2)  Multiplied  the  white  noise  by  the  appropriate  image 
data  standard  deviation  when  necessary. 

3)  Converted  that  matrix  into  a  displayable  image  tile 
using  subroutines  SCALE  and  INTBYTE,  when  necessary. 

4)  Read  filter  parameters  from  an  input  file  into  an 
array . 

h)  Implemented  the  equation: 


1-1  , -1 

y  (  n  ,  m )  -  -  a  j  j  y  (  n  -  i  ,  m  -  i  )  *  w  !  n  ,  m  j  (  A  .  \  i 

l  0  ’ 

(  i  ,  i  )  /'> 

using  the  white  noise  array  and  the  filter  eootticient 
array  as  inputs. 


6)  Converted  the  array  result  from  4)  into  a  displayable 
image  file  using  subroutines  SCALE  and  INTBYTE. 

7)  Used  subroutine  SUBINTFILE  to  create  image  data  files 
from  filter  results  for  further  processing. 

8)  Used  subroutine  NONC  to  apply  the  summation  filter  to 
image  data  when  necessary. 


NONCAUSAL  (VAX/VMS  FORTRAN  1 


This  program  did  the  following: 


1)  Read  filter  coefficient  values  into  an  array. 

2)  Read  image  data  from  an  input  image  file,  converted  it 
to  integer  values  using  subroutine  BYTEINT,  calculated 
the  mean  from  the  data,  and  placed  the  data  into  a 
real  array. 


3)  Implemented  the  equation: 


y ( n , m) 


L  L 
i  =  -  L  ]  =  -  L 


a j  j  x ( n- i , m- j ) 


(A. 2) 


using  the  image  data  array  and  the  filter  coefficient 
array . 

4)  Called  the  subroutine  NONC  to  implement  the  equation 
in  Step  3)  a  second  time,  when  necessary. 

5)  Converted  the  result  of  Step  3)  to  a  displayable  image 
file  using  subroutine  SCALE  and  INTBYTE. 


This  program  performs  the  same  basic  functions  as 


NONCAUSAL,  without  having  the  capability  ot  calling 


subroutine  NONC.  It  was  used  for  convenience  in  convolving 
certain  filter  structures  with  certain  test  images  directly. 
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4.  MAKFIL*  f VAX/VMS  FORTRAN) 

This  family  of  programs  was  used  to  create  various 
autoregressive  and  FIR  filter  coefficient  files,  using 
source  data  manually  entered  into  the  program. 

5.  SPECOR2  C IBM  SYSTEM/370  3033  VS  FORTRAN  1.4.1) 

This  program  implemented  the  equations  derived  in 

Appendix  B  and  created  data  files  used  in  developing  the 
corresponding  graphs. 

6.  SPECOR3  f  IBM  SYSTEM/370  3033  VS  FORTRAN  1.4.1) 

This  program  implemented  the  equations  derived  in 

Appendix  D  and  created  data  files  used  in  developing  the 
corresponding  graphs. 

7.  SPECOR3A  f IBM  SYSTEM/370  3033  VS  FORTRAN  1.4.1) 

This  program  implemented  the  equations  derived  in 

Appendix  F  and  created  data  files  used  in  developing 
corresponding  graphs. 

8 .  VARIMGS  (VAX/ VMS  FORTRAN) 

This  program  was  used  to  display  image  data  files  on 
the  COMTAL  (not  written  by  author) . 

9.  PIECE  (VAX/VMS  FORTRAN) 

This  program  was  used  to  make  128  ■  128  image  data 
files  from  larger  image  data  files. 

10.  INTFILE  (VAX/VMS  FORTRAN) 

This  program  created  appropriately  formatted  integer 


files  from  input  image  data  for  further  processing. 


11.  TRANS  (VAX/VMS  FORTRAN) 

This  program  changed  the  format  of  filter  coeffi¬ 
cient  data  files  into  a  form  readable  by  the  image 
processing  programs. 

12 .  NSHP  (VAX/VMS  FORTRAN) 

This  program  was  used  to  convert  quarter-plane  auto¬ 
regressive  filter  coefficient  data  to  non-symmetric  half¬ 
plane  autoregressive  filter  coefficient  data  based  on  the 
transformation  outlined  in  Chapter  II,  Section  D. 

B.  SUBROUTINES 

1.  PGAUSS  (VAX/VMS  FORTRAN) 

This  subroutine,  written  by  C.W.  Therrien,  was  used 
to  generate  zero  mean,  unit  variance  white  noise  using  RAN 
(a  random  number  generator  function)  SQRT,  COS,  and  SIN 
FORTRAN  functions. 

2 .  SCALE  (VAX/VMS  FORTRAN) 

This  subroutine  takes  an  image  data  array  and 
converts  it  to  an  integer  array  with  values  between  an  input 
maximum  (MAX)  and  minimum  (MIN)  using  the  following  scaling 
formula: 


Ki,j)  =  ,(A(i^±-LOW)  x  (MAX-_MIN)  + 

HIGH-LOW  v  ' 

A(i,j)  is  the  input  image  data  array,  I  ( i ,  j )  is  the  output 
integer  array,  and  HIGH  and  LOW  are  the  high  and  low  values 
of  A(i,j),  respectively  (calculated  in  this  subroutine). 


This  is  done  to  provide  appropriate  values  for  image  files 
that  will  be  displayed  on  the  COMTAL  Vision  One/20,  since 
the  gray  scale  intensity  level  of  each  pixel  is  represented 
by  an  8-bit  word.  So  values  possible  (in  base  10)  range 
from  0  (darkest) ,  to  255  (brightest) . 

3 .  INTBYTE  AND  BYTEINT  (VAX/VMS  FORTRAN^ 

These  subroutines  are  necessary  since  data  in  an 
image  file  are  stored  in  two's  complement  form.  The  related 
variable  type  in  FORTRAN  for  these  values  is  BYTE.  To 
process  image  data  using  FORTRAN  implementation  of  the 
appropriate  formulas,  these  byte  values  must  be  converted  to 
integer  (and  eventually  real  using  the  FLOAT  function)  form. 
Results  of  image  processing  formulas  in  real  form  must  be 
converted  to  integer  (using  the  INT  function)  and  then  byte 
form  to  be  placed  in  image  data  files.  INTBYTE  converts 
integer  type  variables  to  byte  type  variables  using  the 
following  criterion  (I  is  an  integer  and  B  is  a  byte) : 

If  I  <  127  and  I  >  0  then  B  =  I 
If  I  >  127  and  I  <  255  then  B  =  1-256 
BYTEINT  converts  byte  type  variables  to  integer  type 
variables  using  the  following  criterion: 

If  B  >  -128  and  B  <  0  then  I  =  B+256 
If  B  >  0  and  B  <  127  then  I  =  B 

4.  SUBINTFILE  (VAX/VMS  FORTRAN) 


This  subroutine  performed  the  same  function  as 
INTFILE,  but  could  be  called  by  a  program  to  operate  on 


processed  image  data  arrays,  rather  than  just  image  file 
data  inputs. 


5.  NOMC  ( VAX/VMX  FORTRAN) 

This  subroutine  performs  essentially  the  same 
functions  as  NONCAUSAL,  except  that  it  can  be  called  by  a 
program  to  operate  on  an  image  data  array. 

C.  APL  FUNCTIONS 

The  APL  systems  on  the  IBM  System/370  3033  and  VAX/UNIX 
were  used  for  matrix  manipulations  and  operations,  for 
graphing  filter  structures  and  convolution  results,  and  for 
calculating  autoregressive  filter  coefficients  from  image 
data.  All  APL  functions  except  MAKMAT  were  written  by  C.W. 
Therrien. 

1.  MAKMAT  (IBM) 

This  function  was  used  to  create  the  large  coeffi¬ 
cient  matrices  (A)  used  in  calculating  the  FIR  filter 
coefficients  as  outlined  in  Chapter  III. 

2  -  CC2  fIBMl 

This  function  was  used  to  circularly  convolve  the 
Laplacian  FIR  filter  and  its  various  inverses.  Appropriate 
zero-padding  of  these  filters  makes  the  resulting  circular 
convolution  equivalent  to  linear  convolution  [Ref.  2:pp.  VO- 
72],  which  was  the  desired  operation. 

3 .  GFTDATA  (VAX/UNIX^ 

This  function  is  used  to  transfer  image  data  files 
' •  ' : n  a  UNIX  subdirectory  to  an  APL  workspace. 
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4.  PUTDATA  (VAX/UNIX1 

This  function  is  used  to  transfer  filter  coefficient 
data  files  from  an  APL  workspace  to  a  UNIX  subdirectory. 

5.  MEAN  (VAX/UNIX^ 

This  function  is  used  to  calculate  the  mean  of  an 
image  data  file  for  use  in  the  APL  function  COVF. 

6.  COVF  (VAX/UNIX) 

This  function  is  used  to  calculate  terms  in  the  2-D 
covariance  function  for  use  in  the  APL  function  CORR. 

7.  CORR  (VAX/UNIX^ 

This  function  is  used  to  estimate  the  2-D  covariance 
function  of  the  image  data. 

8.  MVLEV  (VAX/UNIX) 

This  function  is  used  in  APL  function  FF2DLEV  to 
calculate  necessary  parameters  for  the  2-D  Levinson 
recursion  from  the  covariance  function  of  the  image  data. 

9.  FF2DLEV  fVAX/UNIX^ 

This  function  performs  the  2-D  Levinson  recursion  to 
solve  for  the  filter  coefficient  vector. 
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AUTOREG 


C  THIS  PROGRAM  GENERATES  A ’l  IMAGE  TEtTlRE  MS  I  MG  sHITE  “OISE  AS  A  v  INPUT 
C  TO  AN  AUTOREGRESSIVE  FILTER  *MnSE  oaraom£TERS  ARE  Ob  T  A  [  NF  D  FRO-  TH£ 

c  fi l£  fil:oef.  subroutine  pgauss  is  used  to  generate  the  input  ««ite 

c  NOISE  ano  SUBROUTINES  SCALE  AND  I  n  t  B  r  T  E  ARE  USE  0  TO  PREPARE  IMAGE 

C  0  A  T  A  ARRAYS  FOR  0 1 SPL  A  Y  .  SUBROUTINE  SUBInTsIlE  IS  JSEO  IF  AN 

c  INTEGER  FILE  RESULT  (S  DESIRED,  and  SUBROUTINE  NONE  IS  USED  IF  A 

e  L*3LAC I  An  INVERSE  FILTERING  STEP  15  NEEDED, 
c 

C  OEFInE  variables 

byte  a(  3: 127) ,bi»(0: 1 27,0: 127 ) 

integer  n,seei«riite>:site>o*l,a*tl,i,j,iol,roa,col,(nteq(0:)27,0:i 

:  2  7 ) 

c«a 1 »  9  yell, v»T2>yn(0:I?7, 0:127), jsU'»,i*(o:|27,Ci:i27),coet(0:b,0:<> 

!),Y8r,ms«,yiin,tnOut(0:l?7, 0:127) 

C  OPEN  FILES 

O3en(jnitsl,ne»es'Cretnii«-»o.0atalyn.aet‘,t»pep‘'<e-,»ecCe*s*‘bireet 

I  ’  .recorisi z«=J2,ma«rec=l2A> 

03ti(jnit:i,n)»*:‘  I  r»t»»Ann, i»l  TioalRa.det ',tv3f:'nt«,,»ccess:'t)i 
!ftCt,<r*C0Miilt:5?,*a>r*::|2o) 
ooeo(5»Tile='  If at lit Jl  *c30u3.3»t ',jtAtus:'oia') 
c  OEFInE  parameters 
*eed=t2jaSb7 
r$  '  Tf: 1 2  7 
c  s  i  te  s I  27 
o  «  I  s  3 
3«1:5 

t>a«  s2?R.  3d0 
■»  i  ns2  7  .  Ort it 

c  CREATE  a-UrE  NOISE  array 

03  10  isO.fjite 

do  20  i:0,ct«tfl 
j  R  1  s  i  ♦  1 

cell  oqeue s ( seed, y a  1 1 , y a  1 ? ) 
en(i»i)=yall 
no ( i , i 3  I ) s va I  2 
20  Conti nue 
10  Cont  tiu; 

c  SCALE  array  a  N  A  \  0  CONVERT  T  u  BYTE  Form 
C  call  s c a  I e ( ■ O , i n t e Of *a « , n < n ) 

c  call  i  n  t  ov  t  e  I  i  *»t  eR,  b  i  * ' 

C  aR  I  T  E  THE  .<m<TE  NOISE  IMAGE  array  TO  A  FILE 
c  do  JO  i  - 0 ,  f  s i * e 

c  oo  R 0  i=0,;JUt 

C  s(  I  I  sot  n( l , j  ) 

C  00  conti nje 

C  yf  1 1  f(  I  ‘  i  t  !  )  (a(n),n:0,fsucl 

C  JO  co-it’nje 

c  REaD  filter  paja-etERS  I'.'TD  an  array 
do  50  isil,9*l 
do  SO  j  s 0 , o *  1 

read  I  J  .  55  1  : or  * (  i  ,  i  ) 

55  l:r«at| o?0 ,12) 

pO  C  0"t  t  1\J4f 

50  Conttnu. 

C  MULTIPLY  a  -t  I  T  £  noise  3  Y  REAL  I  'AGE  STANDARD  GEvItTION 

yarssorf(coeTlO,0)) 

3  o  TO  isu.rjiae 
do  90  jso.csue 

wol  i  ,  j  )Syl(  l  ,  i  )  • V  af 

00  c  ont  t  oue 
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•  1 
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AUTOREG  (CONT.) 


70  c  ant  1 1U* 

AP»LY  NOISE  TO  The  a  mSEGFESSm  filter 

aj  110  nsn  irsin 
ao  120  »>=0*;sirp 
osu»-0 .000 
33  1  30  »  =0  .  on  I 

O'  1  u "  i  A  1 '  .  3 -  I 

» < ( ( l ,ri. o ) . ano . ( I .eQ. 0 ) )  33  is  ion 

r  3  •  -  n  -  i 

C  3  1  =■»-) 

tt( (roe. It. 0). or. (col. It  .0)1  33  (3  1  R  0 

3Su"-s(C3tl  1  1  ,  j  I  *  i  *(  ro.icol  1  I  ♦  n»un 

1 00  C3"t  I  nut 

1  JO  C3"l  i nut 

i  e  (  o  t  e)r(-1.0o0#nsu*)tenln,e) 

1 ?0  coni i nue 
110  Conti  nut 

filter  the  i mage  data  array  jsirg  the  l*°lacta»j  inverse  filter 

Call  nonC(i'<i»3ut ) 

SCALE  THE  RESULTING  IMAGE  array  AND  CONVERT  TO  BYTE  FORM 
CaM  sctltlmoul  i  intt3itti(»m) 

Call  intOvtf{in(C3ini'l 
*  R 1 T  £  THE  GENERATED  IMAGE  INTO  A  FILE 
Oo  1 53  i =0  »  rs i re 
Os  160  i=0,:sue 
®  { i  >  — o  »"><•»#»> 

160  continue 

«r i t e ( 2  * »  ♦  I  )  ( a (n ) , o  =  0. rs i re  ) 

150  Cjntinu; 

CLOSE  FKES 

e I ose  <  un i t  s 1 ) 

Close(unit=2) 
c  I  o  n  e  3  ) 
s  t  03 
end 


NONCAUSAL 


THIS  °  R  3 .1 R  s  4  GENERATFS  an  Image  TErT/Rf  USING  a  n0nC1US»l  fir  cJlTER 
n-OSE  RARl-ETERS  are  O'llAltEO  from  a  Data  aIlE.  The  FIlTE0  IS 
: °3L I E  3  Tu  A '  i  image.  SJHRO  JUNES  SC*LE,  INTByTE,  and  ByTEI.T  ARE 
jSEO  TO  3R£Rape  I maGE  data  arrays  c03  OISPlat.  Subroutine  NO’.C 
IS  USED  IF  At  intermediate  filtering  STE3  IS  DESIRED. 

DEFINE  variables 

Oyf •  a  (  0 :  1  2  7  ) ,  b  i  m  (  0  :  I  2  7 , 0  :  1 2  7  ) 

i iteRtr  n,rsire.csire.i.j.ro»,col.inte3C0:127.0:127),tsire.n(3tre, 

•  inde«l.m3e«2 

r.al  *1  nsu'.  i  *  (  n :  1  2  7  »  0  :  l27),corf(  -10:  10,-1  0:  101,ysu»,.ean,3j*3(0:  1 
1 27,0:127), hiqn,  I  o  «  »  t  r  y  (  -  I  :  1  2  B  ,  -  1  :l?<),«ti,«iri 
0PEN  FILES 

03tn(jni|:|,ni«.:'  (rltlmnn,  i  .1  oar«.  3at  '  i  t*3t:'ol!l'  ,»CCtSJ:  '3'  n 
! e  c  t  '  , rtcorAAt  tt:17 ,«»<r*c :1 J* I 

3Den(unit=2,naee='  Irjtnjnn.  ml  ari  »6 A . 3a  t  '  ,tvoc:'neu'  ,access  =  ‘3t  re 

!  r  t  ',rt;Or3ti  l»:!?i«t<r*c:l?fl| 

OoenfJ,  #i  It:'  Iratneann.  Ht  tl  "c  I  i  Icola.Rat  *  ,tUtu»:'ol  3  '  ) 


t  . ‘  A.«  4.1 


I  I'l  4*4  I'aI*!  II4  ai.  |>.  iVlIa'tl*1 


,«ivv  V  HTnm  X_’r'lJrX’rV'VTY'Y  riWj  v-v  *-„•  r.  rr 


NONCAUSAL  (CONT. ) 


DEr I n£  parameters 
1  ndei I s - 1 
1  nde «  2 =  I  2 b 
rj i ie- 1 > 7 
cs> ?e=l ?7 

«tsir*=-l*fsi*e 
ha«  s2  55 . OaO 
n i nrO . 0  3O 

REAO  FILTER  PARAMETERS  INTO  AM  ARRAY 

dn  10  i  zmiti  it,  Isi  ze 

da  ?0  j :»( si ie» Mi It 
r ea  3 ( 3 , 25 )  cofl  (  i  •  i  ) 

25  t jr »at ( d20 .12) 

20  c  an  c i nue 
1 0  c  a" r 1 no» 

READ  IMAGE  T  3  BE  FILTERED  INTO  an  ARRAY  AND  CONVERT  TO  INTEGER 
da  30  iiO.rsize 

r;ad(l'i»l)  ( » < n ) , n: 0 , 1 2 7  ) 
da  'JO  j  s  0 ,  e  s  »  z  e 
0 i  m ( i , j  )  =  *  C  j ) 

«0  C3»t i iu( 

30  c  jn t  i nu  » 

Cilt  av t e i n t ( 3 i m, i n t e a ) 

convert  rnr  integer  array  into  a  real  array  and  compute  the  mean 

nsun=3  *  OdO 
F  ;  cj«  =  -l 0OOO.O3O 
1 3-: 1 0000. O-fO 
d  3  SO  i:0,r5ife 

da  SO  ji  3  ,  c  <s  1  J  ? 

ini  j  ):(  1<n(  (  i  n  t  ;  a  (  i  ,  i  1) 

i  M  i  » (  i  ,  i  ) .  at.  •>,.»■>)  1  ,  j  I 

'  *  (  t  m  (  1  .  i  1  .  1  t  .  1  -■ »  :  1  3 .  -  1  - 1  1  ,  j  1 

so  C  3nr  1  30; 

SO  c  anf 1  no* 

n;anSnSjm/(12a.0*l2H.3) 

«ri  t;(  *,S5)'if  jn.hi  3(1,  lo. 

55  ijrjaf C  a  r i 3 i n  j I  i*aae  Jjt  1-  m;an  = • , d 1 ? . 5 , •  h » 30  =  ’ » dl 2 . S . '  13-=': 

*  1  ’ .  a  ) 

suAiRacr  me  mfai  cru  t-ie  au; 

da  TO  i  =0,  1  21 
da  SO  j  =  J , 1 ?  7 

t  * ( j , j ) si * ( t .  i ) •«*11 
BO  c  an  r i no; 

70  c  an  t 1  no  • 

r  I  _  T  E R  ME  MAGE  43R.W 
33  R0  nxind;«l.ind;«2 
da  100  *= 1 nle« l « inn;* R 
3su  *  =  0 . 030 

3a  110  i  =  *  M  i  y;  .  M  i  *e 
da  120  i=»<iiie,fiir» 

r  a  ^  -  n  -  i 

Cal  *•»“  j 

iUfrjK.  If.  Ol.op.fc3l.lt.  0  1)  aa  I3  12  0 
iM(ra*. at. 127). or. (cal. or. 1271)  33  to  I  2  0 
d5ir"=f  i  "Ira.icol  )  •  c  o;  <  <  1  f  i  )  1  ♦  3  So* 

20  c  an  r 1  no; 

10  canrinoe 

00  c  an (  i nu; 

R0  cjnr i no ; 

CALL  SjBROjTINE  T 0  1'IVERSE  FILTER  FlLlEREO  MAGE  ARRAY 
3*1!  Annc(trviir) 


'.'  W.V.V  VA' V.V  .V. ‘.'.N  '.'.V  .V  ■ 


V  *-  *.  s  %  *.  - 


NONCAUSAL  (CONT.) 

c  400  Ih£  -e»n  of  Ih£  10  put  IM4GE  Tf)  Th£  fItr£?EO  RESULT 

33  1  a  S  i=0,riif. 

3:  1  S  3  i  -0 , r * «  f • 

i  *  (  i  #  i  ):t  ry ( i f  i  )♦»•»« 
l  SS  c :n r i iue 
U5  Cantnu? 

C  C^rArt  A  ^  fHAGf  )  4  T  A  FILE 

c  *  I  1  suo'^t  H  le(  m) 

c  SCALE  Trt£  *ESJlT  I?4S  l  ^  a  G  E  4  P  ?  A  Y  a  NO  C0vv£P  T  T  3  By  TE  fQpm 
call  sc  i'e<  mi*  i'iteg#"i*i,iMn) 
call  ^t&»te(intnoi  n) 

C  *PfTE  TH£  GEVE3ATE0  I H  AG£  INTO  A  FILE 

Oo  ISO  i s  0  9 r* i t * 

Id  160  i  -  o  #  c  *  »  *  e 
a( 

1  60  c  an  r  i  iue 

writer?'**!)  (a('i)#o  =  .l,rsize) 

ISO  c an t i 3g • 
c  CLOSE 

cioaelunir-)  ) 
c  l  sse  (  u*i  i  t  ) 
c » os- r b) 
s  r  03 
eii 


CONV 

c  r-qs  PPOGPS'*  G£*l£3ir.£S  «  o  1  y  i  GE  TErfjPE  USING  '  ':0OC4'J?4L  P I  a  fILT£a 
S  1*151  3»»*,'£TEas  3-i  0?H|vr>  Fq?*  1  04Tc  rlL£.  Th£  r[Lrra  Is 

c  i0:,L  t  £  3  TO  4-i  I  - 1 GE  .  S  JSPO  JT  [.0£5  S:4lE,  I'.  I  3  t  I  f  ,  4ou  3 1 1  £  I  N  I  *P£ 

;  j  3  7  D  10  » a  £  p  4  a  r  !-« a  G  £  o  4  r  4  aaaars  7  0  p  r>  [  $  p  ■_  »  » . 

c 

C  0E71'-E  v*.0I4aL£3 

oyt.  j ( 0:i?7),o>-(0:l?7, 0:1^7) 

'  n t  eq-r  i,  ps,  „,csi  jt  ,  j  ,  j  ,  r(,<l0l  ,  , ->t  •1(0:  I  ?7,  3  :  I  ?  I )  ,  t  s  , .f  s  ,  je 
r.ai'f  Hull  n>C0:i2?,0:l<?7),C3ef<-M:ll.  -11:11),  isu*.«ta'>,  3^3(0:  I 
:  2  7 ,  C  :  1  2  7 ) 
c  C3£N  FILES 

33fT(jnit:|,4,«.:'  I  ra ! n.sin , 1 ,|  ,r i»5J . , vrf : ' 3 |  V  , acc PS5: ' 1, r . 

:  c  f  l,T-c3risir*:ii,*a<',e::l?,J 

33ei(jnit:J,4)«»:’  Irathulin. i»] jM»s9.]9t  1 »t yoes'ie. ’ ,acc«ss:‘3i r* 

:ct',rrciT3SiT?:5E,Taif(;:li«) 

33*M  j.  *  1  lf='  IralUMm.JJtjICjnvJ]  .  1a  t  '  ,St  atujT'olo’  ) 

c  0£  7  I  E  PiU4“£rEPS 
r  1  1  z  ?  :  1  ?  7 
c  s  1  1  e  =  I  ?  7 
*  $  1 tt-  1  I 

c  PE  4  0  FRTEP  P4P4afTEPS  I'lTO  4N  40PSY 
03  10  i:'»si?»,(jijf 

33  ?0 

r • »3<  3 , 25  t  : se  H 1 , j  ) 

?S  1  jr  »at  (  -i£  '  .  1  £  ) 

20  Cdht 1 ige 

10  c:nn  lua 

c  P£40  I'i'.E  F3  PE  F  I  L  T  £  3£  0  I*JtJ  4  •;  [P)K  too  GDVVFPT  T  f!  [  0  T  f  G  £  P 

33  30  1  * 0 1  r a  1  / a 

rtiK  1  ’  1  *1  )  U(  3)  ,3:3,  I  ?/) 

33  uO  i:“,:su; 
btK(i,j)4a(i) 
a  0  continue 
30  C33| llu* 

Call  })init(3l<i  UM(;| 


87 


■i 

''A 


•\s 

•Vi 

>1 


V'.'.'-’y 

^  1  A  \  ^  *  *  n  M  M  ^ 


A.  - 


•  ■***,  "  *  *  n  *  «l"  ■  "  •  ' 


CONV  (CONT.) 


c  COmv-PT  Tnr  Ala l i  J  N  r  D  *  P  £  A  l  a o q a  y  4^3  C  U  *  p  U  T £ 

^3  SO  3 :0, rj 1 (f 
03  30  l=[l,csl  If 

•'ll'  i  )S'lo»(l  iMfOll,  i  )) 

"  S  u  •»  =  I  n  (  l  ,  1  )  ♦  n  a 

30  C3"( I 3Uf 

SO  Cinmu. 

lost  (  (rji/fl  >  •  (  c  S  i  *e*  I  )  ) 

»  n » tf ( •  ,  S  3  )  »e  an 
SS  131-MI  (  l20.  12) 

c  S'J»t3a;r  Ij£  m  rwrn  N£  I  <«r,£  *3Bar 
33  70 

(13  « 0  i:u,cit;e 
t  *  (  i  #  i  1 

OO  C3flI ' 1U^ 

70  C33( I 3U. 

0  ~  I L  f  E*  >-|£  a -?ai  r 

33  30  n : 0 ,  r  % i  / e 
33  10  3  "  a  'I  *  ;  s  i  Z  e 
3  S  u  -  -  0  .  I  cl  .J 

33  11°  >=*lj!7t,l$ije 

33  1 ?  0  i:>ttl((.ll1|« 

r  3  a  z  n  -  I 
C  3  I  =*>-  i 

» f < ( r 3M . I  t  . 0 )  . or .  ( e 3  I  . 1 f . 0 ) )  33  t  3  120 

i*((r3«.3t.l27).Or.(C3l.3t.l?7))  30  to  1(?0 

d5U':l  I  "(ro-.col  l*C3«l  (  i  ,  i  )  )MSu* 

120  Continue 

110  C3"rnut 

dji3(i,*l:asu> 

100  C  3  n  t l  3U» 

3  0  Confiiu* 

C  SCALE  rn£  JESJlTIM  I  <AGE  tRItr  A\1  CO'jvFBT  T  0  Stic  POSI 
ca''  sc  jl*(iu3o, i 3tea) 

Call  mto»t  t  (  inteson) 
c  ME  GE'iEBAfEO  HAGE  INTO  «  FILE 

33  ISO  i : 0 . r a i z? 

3  3  130  idiCSWt 
a(i)s3i-"(i,i) 

IbO  C  3nr  l  nit 

-rite(2‘i*l)  ( a ( n ) , ns  0 , r J i z;  ) 

ISO  C  3n  r  i lu; 
e  CLOSE  -  I _ £ 5 

Cloaelunitsl  ) 

:  i  ostium  = 
c l oae ( J ) 

St03 

end 


E  «E  AN 


■MMWW  vwmuuww 
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tt  rr  t  .f  \mw\t*  »."  ig*  vnr«? 


-•*YWYm'WX7rlKjrJir<  W/TJV 


SPEC0R2 


C  THIS  PROGRAM  SOLVES  EOLATIONS  FOR  SPECTRAL  CONTENT  AMO  CORRELATION 

c  IN  One  direction  of  a  given  autoregressive  image  model,  it  mrites 

C  THESE  RESULTS  TO  OEVICES  5  AMD  «  RESPECTIVELY. 

C  DEFINE  VARIABLES 
INTEGER  t.X 

REAL *8  Ft, ALPHA. THETA. M. ALPHAS. SXU.RXK<-4», AT ) ,2. ALPHA* 

C  DEFINE  PARAMETERS 
PI  *3 . 141S92.54 
ALPHA. 0. 

ALPHAS* ALPHA*. 2 

C  START  X  AXIS  LOOP  AND  DEFINE  X  VALUES 
DO  10  1*0,  ♦» 

N*t-l.O*PI)»t2.o»PI«(FLOATCI)/»».o)) 

C  DO  SPECTRAL  ANALYSIS 

SXH* 1.0/I1.0-C2. 0* ALPHA* COS IN  I ) ‘ALPHAS) 

HRITEC.ISIM.SXH 
IS  FORNATIF1 0 .S.1X.P10.S) 

10  CONTINUE 

C  DC  CORRELATION  FUNCTION  SOLUTION 
DO  23  K.0.41 

ALPHAK. ALPHA. *K 
RXX(X).(ALPHAK/tl .0- ALPHAS)) 

I  —  1"X 

RXKCl )*RXK(K) 

23  CONTINUE 

CO  30  K—49.49 
7*FLOAT(X) 

hrite(4.;s)2.rxk(k) 

25  FORMATIFIO.O. IX .Flo . 15 ) 

30  CONTINUE 
STOP 
END 
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SPEC0R3 


c  THIS  PROGRAM  SOLVES  EQUATIONS  POP  SPECTRAL  CONTENT  AND  CORRELATION 
C  IN  ONE  DIRECTION  OR  A  GIVEN  AUTOREGRESSIVE  IMAGE  MODEL.  IT  **ITfS 
C  THESE  RESETS  TO  DEVICES  J  AMD  A  RESPECTIVELY. 

C  DEFINE  VARIABLES 
INTEGER  i.J.M 

REALM  PI  .ALPHA.  THETA.  ALPHA  J.AL3PA.  ALPHAS.  ALPHA*.  TMOTN.COSCT.M.  A.  B 
■  .CQSCH.SXtal  .  SXM.*.  ALPHA*.  *THET  A,  ALPM2*.  KP2TMP  ,KTMP,  INTERN.  RX*(-*f: 
•4f >.Z 

C  DEFINE  PAR  ALTERS 
PI-3.141512454 
ALPHA«0.f 

THCTA-0.0«(PI/12.0) 

IF(TMETA.EQ.O.O)  tweta«o. 00000 i 

ALPHAS -ALPHA* *S 
AL3P A* ALPHA J • ALPHA 
ALPHAS* ALPHA* *2 

ALPHA** ALPHA* *A 

CCS2T.COSC2 .0* THETA) 

TWOTH*2 .0-TMETA 

C  START  X  AXIS  LOOP  AMD  DEFINE  X  VALUES 
C  DO  10 

C  W*(-l .0-PI )• (2.0-PI ■ (FLOAT (I J/JT.0)) 

C  DO  SPECTRAL  ANALYSIS 
c  a*theta-h 

c  B*THETA*« 

c  cos2h*cosc:.o*w) 

C  SXH1 -ALPHAS- (COS( A  )*AlJPA)*(ALPHAS»CCS2T )- (CCSIB )*AL5PA )♦ (ALPHAS 

C  ■•CCS2N) 

C  SXH*1.0/(i.O*ALPMA**(2.0"SXWl )) 

C  WRITE (3. 15 JW.SXM 

C  15  FORMAT (F10.5.1X.F10.S) 

C  10  CONTINUE  • 

C  DO  CORREwATION  FUNCTION  SOLUTION 
DO  JO  J*0.4« 

X*PL0AT( J  ) 

ALPHAX*ALPHA**J 

XTHETA*X*THETA 

KP2Tmr.( (2.0*X)«THETA)-PI 

KTHP-XTHCTA-PI 

INTERN* (COS (KP2TMP  )- (ALPHAS ■ CCS (XTMP ) ) )/(l . 0* ALPHA*- ( 2 . 0* ALPHA 
■S*COSCTWOTH) J ) 

RXKIJ  I*  ( ALPHA*/ (2 . 0*CSIN( THETA  )**2  )  ) )  *  (  (COS  ( KTHET A  )/( 1 . 0-ALPHA 

•  s) )  • intern ) 

GO  TO  3* 

C  35  RX*(J)*(ALPHA**(2.0**n/(2.0*((1.0-AL3HAS)**2)) 

34  M—  |*J 

RXK(M)*RKK( J  I 
30  CONTINUE 

00  40  j*-4«.4l 
2  -  FLOAT ( J ] 

ur:Tc(4.*5)2.RX*(J) 

*5  FORMAT (FlC.lilX.FJ0.15) 

*0  continue 


SPEC0R3A 


C  THIS  FROGRAN  SOLVES  EQUATIONS  FOR  SFCCTBAL  CONTENT  AMD  CORRELATION 
C  IN  ONE  DIRECTION  OF  a  GIVEN  AUTOREGRESSIVE  IMAGE  K3DCL .  IT  MIITCS 
C  THESE  RESULTS  TO  DEVICES  3  AMD  4  RESPECTIVELY. 

C  define  variables 

INTEGER  I.J.H 

REAL**  Rt.A.B.AS.AMBS.AS.tS.AlCI .AA2.fR] . BK2  . INTERN. RxR< -A  9 r A  9  ) . 2. A 
•  R  .  RF2  .  RF  1  .  OF  AS  •  OKAS 
C  DCF  INC  FARAHCTERS 
FI*S. 141*92454 
A*-0  .  9 
B— 0  *9999 
A*«A"B 

AKBS • ( A— B  )  **2 

AS*A>*2 

*$•*••2 

C  START  X  AXIS  LOOF  AMD  DEFINE  X  VALUES 
DO  10  1*0.99 

W»(-l  .  0*F  I  )*(2.0*FI* (FLOAT ( I 3/99.0) 3 
C  DO  spectral  ANALYSIS 
C0S:n*CCS(2.0*W) 

SXWl»l .0-(2.0«(A*B*(AS*83*CA**S) )*COS(W) >• (2 . 0*A*B*COS2M) 

SXN»  1 . 0  /  ( SXN  l  •AS*(2.0«A*B)*SS*(AS«BSn 
NRITEC3. 15  3N.SKN 
15  FCRHAT(FI0.5 . 1X.F10.S3 
13  CONTINUE 

:  DO  CORRELATION  FUNCTION  solution 

do  :o  j-o.Af 

IcCA.E0 .33  GO  TO  35 
Aa!«a««C J*l 3 

ax:*a««( j*: 3 

BR 1 »*•• l J* 1 3 
8K2«B**( J*2) 

INTERN* ( AR2/U.0- AS >-(  MAR ]-B)*a*BR:)3/(1.0-AB))MBRO/(1 .O-BS) 

■3 

RXKU  J*(  1 . 0/AMBS  3 •  INTERN 
GO  T0  54 

35  AR*A*«J 

RP2*FLOAT( J*2) 

RFI*FLOAT(J*l ) 

OFAS* 1 . 0*AS 
Ohas*1.0*AS 

RXR ( J ) • ( (K?:«RF1*AR3/(2.0*ONAS3  )-( (FLCAT(J)»RF1*AR)/::.3*0FAS3 

•3 

54  N«— l • J 

RXK IN) « RXR ( J  ) 

50  CONTINUE 

DO  45  J*-49.49 
2 'FLOAT  I J ) 

WRITEU.4S  )2.RXR(J) 


45  FORNATiFlO. I . IX.F30. 15) 
40  CONTINUE 


PIECE 


T-i;S  OJ3SSO’  1  >>•  3  ’  £ :  £  0  ■!  If  *  !^»&r  •'»'*  FI-r 

4  .  ;>  r  *  1 1  0  4  14  i  -j  r  l  4  ■ .  t  *  •  r,  -  --  n.  E 


3£  -  |  -it  *  4-  I  4  sif  s 

o  ,  t  *  ■*  f  D  :  S  1  l  J  .  o ,  •  (  1 :  sit,  o ; s l  I  ) , »o  <  0  :  l  ?  7  * 

Qpc-J  «■  I  L  £  S 

044  4(  Jn  i  t  :  I  ,  n  %  "» •  4  *  (r4t  t*4'>n  »  *  vl  4  4»f  4»Olt  ’  I»ccf4*-  3' 

03f  4  (  j  n  ,  t  Z  2  '  n  %m*  z  '  l?4CDD4nn,.,»l  tfffs44,1|l  *  ,  t  vOf  ;4  «  '  #*CCM*'  3' 

!f:t  'i'e;OMii  ff  :l>il 

HE4.3  I-aJE  77  HE  f  ILT-H?  I'.fJ  »■-  4SS4* 

3d  50  i :0, 1^7 

r*?43(  !*,♦!)  (*(*4),n!i)(137t 

3j  iO  j ;0, 1 ?7 

30  conn  iu? 

30  conii’u? 

,3ii-;  ih;  oesised  !«a3-;  =o»mo-.  into  4  file 
3d  I  S  7  1=5),  IJo 

O  D  167  |  =  6  3  »  l  ^  6 

«  -  I  -6  5 

j)(,  |:3I'(  I.  1) 

1  50  C  D  d  t  ,  *,u? 

«Pl(H/’l*4?)  (  J3(n|  ,  45l,s!l 

ISO  confnu, 

close  -i.es 

c laie(uiir:i  ) 

:  !0,CIU'i!:r’) 

4  r  3  3 
fnl 


INTFILE 

THIS  r-’E'TES  4  1?S<I?H  Q3  5JI51  fjr£j£3  F  l  (.  E  FJO'  ? 

IM3ijr  !'<4GE  cI'.E 


•3  £  -  I  -V  e:  V43M-LES 

o,(;  *  i  .7 : 1  2  7  ) ,  o  • « (  0  : 1  ?  7 , 0  ;  1  ?  7  )' 

int;a(P  i(i,<',«,l,int?a(0:l£7,0; 1 ?7 1 ,4J»>5>D5, l o« 
n,4 I  *  3  »e»n 

: p-w  files 

ooenl  jni  t  =  l  iO)'f  :  'Iratnnjnn.i-linorJ.iat’npofc'olD'.iccpisF'd 

I  irfct  '  , rPC0n7s I  z?  =  57i  ,4,n*c:175| 

ooenlc,  f  i  ;»='  Iran«di4,ij!al  3  ,  -  «  4  •  .  3  »  t  '/SlutuS:'"?.'  ! 

3£43  I-43E  73  H-  f I L  T  E  ^E  3  17  10  4. 

33  so  i:0.1?7 

Pf#l(  1  '  ,  ♦  t  I  (»(-5)»P5  =  C,l?7) 

0 3  7  0  1=0,  I ?  7 

70  C  5"P  5Uf 
50  c  DHt  I  DU? 

C0HVEF7  t  w  r  a  v  T  E  I \ 7 3  4N  !v7fCf4  »»7»7 

Call  Drrfinf  In'  51  tn(?Dl 


vv~-  rj  r-» 


INTFILE  (CONT.) 


lt;n  4N0  L3«  UU £3  OF  Th£  I^a-.E  Jtn  »N0  aRIIE  ID  T£<»«l 

*  i  it =  i 
l  3.=£5S 
15  51  i  =  0  .  I  ?  7 
1  3  5  ?  /  =  0  .  I  ?  7 

>M  i  ,  ;  l.l(  .lo,l  li»:ii(«](i,il 

'  1  (  illlll  i  ,  |  I  .5!  .11  j«  |  M5i:i1l  fl(  i  ,  i  ) 

5<?  CDHI  1<J» 

51  C5"niu< 

■  rlffl  •  ,  Jilhljl,  1.3  m 

55  ♦ jr  »*  t  f  > i o ) 

CDvOjl;  -,£»M  OF  r  N  T  £  J  £  ^  INI  St)5T’4C7  If  FSn-i  Th£  Dili 

U":1 

10  5  5  > I D . 1 ?  7 
io  d  3  |S'i,  l?7 

5U*:SU1*1Atf3(  I  t  i  ) 

5  5  ;  o  "i  r  i  n  j» 

n«a->=<  I  jar  (  *u  «  )  /  (  I  ?1 .  o  •  I  ^  .  0  ) 
wnt;(  a, 

5  j  tor«»t(ilO.?lC.5) 
io  5 i  <  s 1 » 1 2  7 

13  J-s  j  m  ,  I  ?  7 

: o ' t  \ r»ue 

5 1  :sm  nj* 

,  i  r  r  f  l'ii£5E:  i^>>  ivtj  t  i?5<l?f»  cm  f  j  L  c 
,  =  ia-< 
i  =  ia-> 

.r 1 t • ( ?,  5’  1  N  .  I 

53  Hr«M(Ji5) 

30  S7  1=0, 1^7 

33  07  i = 1 » 1 1 ? , 1 3 

•'ttsil  ,  la'll.  i  i«3l>  "M*3l  '  .  i«»>  •  'MMIm  ill  I  .  1 

,iff3(l»i*7),iiri3(i,l»l'l).'iTa5fl,,»ll),lit*3(l,i*17),latf3ll 
'  't  fll  '  I  i  '  I  1),  llt»3(  '  ,  i  all) 

7  i  »da««t ( 1 jii) 

on  : 0 1 t l ng( 

57  ;31( ' 1 jf 

CLOSE  -i.e- 

Clllflui'fFI ) 

: i oseluii i :?) 
s  r  03 

e  oO 


>,j)»'7tfl(l/i7l),>7te3(i,ia?),ti'»e-3.',, 

•ill  «nt*3l  >  ,  i«»l  .  17(13(5  i  i«7  I,  .  I 

17).’3la5fl,l»ll),li(*3fl,i*17),lr[a3li 


i  *5).  ' 

♦  r  \  ,  \ 
,i*15) 


TRANS 


T  H  r  S  0£  4  )5  0471  our  IF  4  FJS^r  IN  F3£r  FOUHT  450  C0NV1: 

17  TO  4  O90G94*  RE4046LE  Fo-"4!. 

1  it  f3?r  i  ,  i 


real  •*.  a  (  I  0  ) 

03Ci(  1  .  M  I  f  =  '  Inimm,  3at  a  I  a  3«a 
03fi(Ji<i  I  e  =  '  |rj|mni,ij(i|  <;i. 

f e»3( 1 >  >)i » i( l( l ),s(2), i(5) , a(4' 
•  ll).i(|7).l(15),i(m,a(l7).KI 

do  10  * :  1  ,  1  6 

mr i t  ? (2 . 20  )  3  (  k  ) 

20  l;a-it (a?0. 12) 

|0  C  3"»  <  INf 
St  03 
nd 


. .  .  - - - »t .1st  * .s»  stuss'ald* ) 

<?  >  t  >  i  e  = '  Irjimni.mii  ':i.it,n'  '.sinusF'ifa1  ) 
li*)iii,5(ll,i(i),j(5),aU!,)l5),al6)ia(7),i(fl),a(d),a(10),j( 
(!?).»( 15),  add),  »(  1  5)  ,a  (  i  6  ) 
i  ■  1  ,  1  0 


$0v 


iV\>\i>/,ivj 
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THIS  GE*lE=sTE3  M  MiGT  T  £  *  T  (Or  |  S I  ^  -  .hit:  NOISE  4S  » '<  INPUT 

TO  *  \  S'JTCIPEG’ESS  I  .'E  F  II  IE  5  |H75E  31H3-.-rr35  40.;  nHT4H.EC  '50-  T-; 
:[LE  F  t  l  :  OE r .  SJo»n  IT  I  \ic  OGl'iSO  IS  JSE''  TO  T«;  If;PJT  ,h]T£ 

-jOISE  »N)  S'J*S0UI  InE3  SC4LE  4N0  I  y  T  B  y  I  e  »»E  iSEO  10  Pocp*S£  I^iGE 
•)4T4  43*47$  F03  DI$?L*T. 


OE-I'JE  v4SJ4?lE$ 

o  *  t  ;  iU:ii7).o.*(0:l?7,0:l?7i 

iittOff  n<Jfej<oiI*<:4ilf>o,'liT*l<'ii>iol>ro»<col,inteo(0:lE71(l:! 
! 27 ) » Horn) 

r.al *6  *all,vsl2...n(0:l27.0:127),  3su'.."<<0:l27,0:i27),coe*(-2:<7.-2 

j  :  o ) 

QP;N  FILES 

03an(jnitol,Pa',?='  Ipatf'i'a'>n,aat4l»n.7at,ttv5e:'''e.,,a1.tess:,airfci 
!  ',reC3r3Si*#452.i**«rer:l2S) 

Osf'l  J''it32,n4»f5'  Irattiuin.iaU)  rofate.lat  *  ,access=‘3i 

0  3«i(  }<  f  i  It:'  lrj|4*ji'p.Tital  fc  Ja.  3a  t  '  i5t9tus:'3l  3'  ) 
o  3  e  i  (  J «  'i  Iffo*  |r  jtnttna,  lat  al  t  s  •  i  a  t  ,iSl»tus:'',t«'  ) 

DEFINE  P404MPTess 
*  *e  3= 1 2  5  JSb  7 

P5'7t:l ?  7 

c s ' re z I  3  7 
o  .  1  =  2 
3.  I  -U 

=• 1 *Ott 

CS;\TE  a  --  I  T  E  NJI3E  4004Y 
03  10  > :0 , r% i ie 

33  ?fl  i=0,CS>7e-l 
i  3  1  :  i  ♦  1 

call  3 oaus s l s tt 3 « V  a  1  1 , n a  I 2  ) 

«3(i»i>=wal| 

.'In  isl  ):vat? 

2  0  e:M  i  iu; 

10  c  3n  t  >  3U  • 

SCiuE  •10  4  no  COnvE’I  T  0  JfT£  p)*'* 

call  scale («ni ' "(to) 

Cal'  mto»ttli'i*3i3i*) 

«0|TE  Trt;  ,»Ilr  'IUISE  I  "4GE  430  4  1  TQ  4  FILE 
03  30  i:0,rsi jt 
33  JO  j  :(',cs>  je 
a(j)s3t'i(*»j) 
a  0  cain  *iu? 

.pilell'itl)  (a(a  )  i':'7i  ts  i  iel 

30  c:'f i in* 

Oil)  FILIE’  =>»  **  p£  TE  3S  l;7)  4N  4SH-.Y 

03  30  I :ss»l |J»I 

33  $  0  f  S  0 i  3  t 1 

p*a j( 3i l  c of t (  •  i  i  ) 

SS  l3p.afI3EC.lE) 

I'll  CJP'MU. 

31  :  3  3  t  •  3  u  • 

1J3L'  .HIT:  3.  _i  ;  5  E  7D  I  hE  4  1 1  0  Jp  G  f  E  ?  $  I  E  rwl:° 

33  117  .  =  Ti,csiJr» 

03  12o  .:0 ,  r s 1 1 e 

OS u.=0.03P 
33  1  JO  I  - -O . 1  1 3- 1 

3  J  I  4  0  j  S  f)  ,  3  3  | 

Mil  1  .eo.'l  .1.(1.  I  i  .pc.  '1  )  33  t  3  1  'JO 

r0«SH-> 

C  3 1  i 

I  |((P3«.  If  .0  1  .OP.lCoi  .  It  .  'I  1  33  I  3  lu0 

3SU~-(C3b‘l li i ) • 1 .( PO.iCOl 5 i'Hj' 

I  J  0  Cl'MIUf 

110  C  3H  f  I  nu» 

■  .  (  -  !  .030.3SU.)  »  h  n  (  n  ,  ■»  ) 

120  CJ'li  VJ( 

110  (3.HHI! 


I  •  I  »  L»  (,» 


.u  ,u-Ji  -«*  t^4 i 


NSHP  (CONT.) 

SCALE  THE  RESULTING  IMAGE  AMD  C3MVE3T  T3  6Yl£  FOP* 

call  sc  j I e  < i * , i nt  eo ) 

Call  (intern) 

tM£  ‘j£»«  rt  1  MAGE  into  a  F  t  L  £ 

3  a  ISO  i'liMiit 
33  160  i - o • Z S I  I e 

a(j)  =3  <■»(>.  j) 
loO  cjnnnuj 

HP! t((?' i *1 )  (j(n),n:0irsile) 

ISO  C3"t 1 nus 
CLOSE  -  I  lE  > 

C lojeluni t  =  l  ) 

C  I  Ose  <  U3«  t  =2 ) 
c I ose ( J l 
C I 0«s( 4) 

St  06 

end 


PGAUSS 


SJAS0‘)TIS£  °0A  jSS(<  .  21  .  ZD 
oE!L»*  a.  3, 2  1,23 
IM1£0£»*j  n 

a  i  SO?T(-?.000  *  4L0G(M»'4f  * ) )  ) 
3  :  S.?851S500  •  -11  NM) 

?i  =  a  •  cosrn 

;2  -  i  *  STUB) 
j£  rj?fj 
£  N  0 


SCALE 

S  jP  r  3  J  t  lit-  3  C  J  I  ?  t  >  n  r  , 

Ti.'S  S  8  1 0  J  1  t '» £  nClLcJ  «i  Hill  M  f  ■ .  T  £  S  c  <  * 
-A  ||<J<  1-0  M  -I  "J  < 

?£-!';£  Y  'a  I  *8L£  S 

>  •>  r  •  3  •  -  iit«3('J:l?7,0:i?7i.i(| 
r  e  a  1  •  3  jnp(":l77(0:!3))-nn,n#l3rf,-ienn,ia< 

=  Elc:)9‘  3  C  4  L  l  3 1 

n  i  am-  1  10/0  .  OiO 
1 3-  =  l 1000.0  to 
33  1 0  i in , I ?7 
33  M  jrO,1?7 

1#(lPn(tij).l».)0«)  I  3«sam(  «  ,  j  1 

t  M  )fpl  i  >  i  )  ,3»,ii  i  jn  i  ,  I  I 

Jj«:»rr(i , j ) 'Sj» 

20  c  3"t i nue 

10  cjntnu* 

I?8.0«1?S.O) 

-ft  '?(  •  .  2S  )  n  l  33,  I 

2S  ^r-at  (  1  ones;* lea  t  -*aq  •  tarn-  m  in:'  ,  i|^ 

!  i  J  i  ?  •  3  ) 

COjfl'ijE  SEAL  I  ‘IS  AND  COni/FII  To  MfEGFi  rj3i' 
33  30  130,127 

33  JO  i  s.1,  1  2  7 

l  3 1  •  1  (  1  ,  j  ):\-it  (  t  (  jrr  (  i  ,  ,  >  -  i  -n  «  ,  •  i  l  -  ,  ,  - 

«4  0  C  ln  t  t 

3U  Cin-ny. 


<  »lLJt  ; 


112.3.  • 


111-,.-. 


AO-A182  144 
UNCLASSIFIED 


IMAGE  TEXTURE  GENERATION  USING  AUTOREGRESSIVE  2/2 

INTEGRATED  MOVING  AVERAGE  (ARIHA)  HODELSCU)  NAVAL 
POSTGRADUATE  SCHOOL  HONTEREV  CA  S  C  RATHMANNER  MAR  87 

F/G  12/9  NL 


SCALE  (CONT.) 


CALCIC*!;  -elGi,  L0*«  A'ly  "£*»  Cf  SC*'.£3  I-*G; 
s  j«* 9 . 0 lO 
*t  3»s-l  3PC0. 0 lO 
I  3«s  1  JOIO.O-iO 
oi  50  i  si)  #  I  27 
3a  *>0  j  s'),  1  ?7 

»nlel3af(i'sre3(i»i)1.3t.,"a*'l  Si^sdnf  (tofejlii  jl) 

1 1  (Hjif  (  mrei(  i  >  i )  )  .  1 1 .  I  lD»s*i9«f(i*»t«a(i»j)) 

60  ejntMu* 

50  ejotHuj 

•»»a-'*^oii/(l^5.0«t?9.0) 

• ri  tel  •  ,  65 )  r*  t  on ,  t  om ,  *«*i 

65  *3r»et  (  *  sc  tlfl  »  <ua*  Oats*  H i a*s • #  3 1 2»5 e  *  lo«z*  ».ll2.5»  1  spans',  jl 

:2.5) 

net  jn 
e-il 


INTBYTE 


S jOno jt in»  i  n  r  o  v  t  •  l  >  "  t  e  3  »  ?  >  *■  1 

3JS  SjdPOjIJNE  TAnES  AN  1<T£GE**  ‘"‘t  C3NVEPTS  IT  I<T0  A 

art;  AJRst 
ErlNjt  \M»IA3LcS 

i  n  t  •  a  *  <■  i,j,i'nt«3(r<:137,!';12T),ni 

Dvtf3tnI0jt?7#9M27) 

cpf'JP'’  CO'ivEPSl ON 

nsO 

33  |0  «  s 0 , 1 2 7 

do  20  j  sO  » l ?7 

iM(inte3f»,j).It.0).3r,(inee3(i»i).3t. 255)  >  «  =  "♦! 
lMintea(i»i).,e.f) 

i  *  t  iKtecl  i  #  j  ) .at  •  255)  3 1  •"  (  i  .  i  )  z  '  2  7 

«M(»->tea(»,i).te.l?7).n''3.<»3te3'>,i).3?.°)>  o « •»  C  »  .  > )  = 

‘  ^t'?<{(2-st-a(«,i).3t,I27).»n'3.(s3teati,i).te.255))  6s»(i,i)' 

! *  nt  »af • , i ) -25a 
0  es"t *  sue 
0  C  3rit  t  nviy 

•  r  i  t  el  *»  IM  *n 

0  IsrsitC  THE  'lU'nSEP  0r  ’Of.'TS  Oil  0C  PANG-.  15  #i5) 

PftuM 


_:s..vt„  In  DT1C  doftK  no< 
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BYTEINT 


ijorsjttn*  ay  f  e «  nt  (  3  i  « .  t  "i  d  ) 

THIS  “ROJRA**  Takes  A  HrlJ  ARRAY  *‘0  EO'‘VFRTS  IT  INTO  AN  INTEGER  ARRAY 
0  £  -  I  'lE  VARIABLES 

tats?**1  i  •*  t  »3<  0  s  1  27 . 0  :  t  2  7  )  .  i .  j  .  n 
Oyt*  3 i » ( 9 : 1 27 . 0 : I  2 7 ) 

PEJFjRR  CONVERSION 
n:0 

3a  10  isO.127 
ds  20  jsO.127 

»f(ft>»»(»,i).lt.-129).or.{3»»(<,i).3t.l27))  asnyl 
ifCj«n(i»j).Tt.*120)  intealii j):-l2S 
if(3in(i,i).3t*l27)  iatpQ(i.j)R|27 

i  ♦  (  (oi  •»(  i  ,  i  )  .qe.-l  25  )  .  and  .  (r»»  *(  i  »  i  )  .  I  t  .0)  )  inte3(i>i): 

I o i n ( » . j ) ♦ 25o 

<M(oin(i,j).3*.0).sn5.(pi»(t,j).le.l27))  inteq(t.j)soi«(i«j) 
20  e  sat  i tut 
10  e  jot i iu) 

•  f  i  r  •  (  •  .  $  0  )  n 

10  Ijr.jtc  HE  VIH9ER  0-  =0IMTS  CjT  0?  RANGE  IS'.IS) 

f»tjf 

e-»3 


SUBINTFILE 


Sj3r3jti"«  Su3>'(<'  If  I") 

TH(5  S.'SROiTINE  CREATE:  a  1  2  *  «  !  2  3  I  N  T  E  G  E 3  Oata  "  h.  c  p  ®  0 «  f  1 1  T  E  A  E  "> 

T  «  4  G  E  0  A  X  a 

0ErlNE  vaRtaa|_£3 

t»>#ri»*»l.(atea<<>tl2T,0;l2T3,*ua 
r  t  *  1  •  H  «*jn.ia(0:i27,;;:|27),ni3a«)3. 

O^EN  PP.ES 

0  3fi(  l(f i  If!'  li'WiJr.ljt  ’  ,StStUS:'"*»'  I 

CO^R jTE  iISH  AMO  10*  VAljES  Oc  ThE  JaASE  OaTA  AND  aRITE  to  TERMINAL 
0  i  3"A-  1 . 033 

I  3«s I  . 0  It 

33  31  I  s  0 « I  2 7 

33  32  jso. 127 

It. l3«)  lo.At«<i,i) 

»  *  ( i  ■»(  i  .  j  )  .  3t  .  ">  i  qa  )  o  >  31s  1  k<  «  .  j  ) 

32  c:nt ' iu* 

31  Cjnt 1 iu; 

«<‘'te*,.33)hi3".l3« 

33  I  jr»At ( ?J2&. 1 2) 

OOhRjt:  5 t  intEOER  array  a»0  subtract  IT  FRny  ThE  OaTa 

>u.=  0 

33  35  i SO. 1 27 
10  “5  i*J, 127 

iu«=iu.t intt j( I , I ) 

<*5  :sK  i  age 
35  c  o a t  »  a  je 

YffSlIjAt  (ju.)/(l2J.  .'>•120.0) 

3u)sja.ae*n 
3'J  t  o  r  a .%  t  (  *  TO. el  0.3) 

33  3h  i so. 127 
33  “•>  i  s.i,  |  27 

i"ie3(i,  ||:nl(3l  i  .  j  )  -  1  a  t  i  ae  »n  ♦  0 . 5  ) 

«•>  ISHtiHge 
3a  :o’t  >  "  if 


COPY  ovail.  H  Drf,  a 

penni*  fuuy  jtBibl.  It;;od, j,;:,;; 


SUBINTFILE 


NONC 

Ijorojt  ini«  no-iC  (  7  r  y  ,  t  »  ) 

c  this  sjs-»ujTr\^  t*<e$  4.-4  i-aje  *0347  avo  ciltE3s  it  «itm  4 
e  mOsCVJSa.  r\Q  *l'_tE5. 
e 

c  DERIVE  V4S  1 43LES 

titeT»i'  o, r*i *e,ss  » /e. i , J »  r©«. eol . ♦*»  fs » »♦ * i te* I oge« l > i oor «2 
r ra  1  •  *  isuM>t*(i:t?7,3:t27),coeTf*IS:iy#~l0:lG),wswii«*«*o,3jm3(>i: 

•128,-l:i29),nigo,lo-.7r/0:i27,0:l2n 
e  OPEN  PILES 

Oseo ( 1  ,  ♦  i  1  es '  (rimum.  Jjulncti  I  eo2 1  .3*  t  *  »st«tu*s'ol  9'  ) 
e  0EC  I  NE  PA3IM- TE^S 
ioot.ls.i 

rj ( it* l ?  1 
en  ie*l  ?7 
♦  s  * ie« 1 1 

«**t*»S-l 

e  PEiO  PJLIE7  °4*4'<e  T£PS  InI7  4.<  43347 
9a  10  i s«* 1 1  it, * s i ** 

03  30 

p*«9< 1 ,25)  eoef ( * .  i ) 

25  »sr.4t(l27.l21 

20  Csntiouf 
|0  Cjnti-iuj 

C  T4<£  T-iE  l  <(P  Jt  4  9ft  4  7  AnO  C3*3<JTE  T  HE  «E*N.  Hl5o,  *N0  L0«  V  4L  JES 
O'lO 

htOOS-IOOOO.OoO 

I 3*s 13090.090 

03  50  1  nnjKi  1 1  inlet? 

03  30  i * » oOt« t • * niei 2 

■7»v.*TPy( i , j )t'iu< 

OioOSTryli.j) 

i ( ( f r» ( i , i ) . I ( . I o.  )  I  jus 7 rv (  i  ,  j  ) 

60  C  30  7  i  107 

50  c  307 1 iui 

«(4isi|  j</)  lg)t  1  (  i,'ie«J»<i|l7<  I  *1  )  *1  lnieiJ'iiOeiMUl 

•rit7(*,55)'74«.ii37i •  O  « 

55  (smtC  toout  i  «*oe  irrit*  *  #  01 2. 5»  *  0  i  o*s  ’  ,  flt  2. 5,  *  lavs', 3! 

.2.51 

98 


‘  -  -  ** 


CO'IVEOT  I  MICE  04  T  4  1MIO  INTEIE3  f  0  PM 
03  tf 3  i s U , I  27 
03  55  I  SO, l?7 

•  07  «3(  i  ,  j  )s,  nt  (  i  v(  i  ,  j  )  .0.5(10) 

55  03074  mi 

OS  C3"t • lut 

*«ITE  IN(E;E3  40347  I *4 T 0  A  130X12?  0414  SHE 
ssI2« 

I  sl?8 

«n  te<  1  ,  I 

JO  * 3P ««7 ( 2i 3 ) 

03  50  isO.127 

03  60  j  so , 112,16 

-r«  r<»n  ,  7  7  1  lo.»3(  i ,  } ).  iotvqfi .  i  .n.  i«teo( , ,  j  72).  «6t*a(« ,  j.  J),i 
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end 
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APPENDIX  B 


DERIVATION  OF  THE  POWER  SPECTRUM  AND  AUTOCORRELATION 
FUNCTION  FOR  THE  TWO  POLE  AUTOREGRESSIVE  MODEL 


Power  Spectrum 


|  H(e3a))  |  2  =  H(e3a))  •  (Hfe^  )  )* 

In  this  case: 

H(ejw)  =  - ±-r-  H(e-jw) 

1  +ae  J 

Calculating  H(e^*)  )H(e"3w) : 


H(eJ“)H(e"3w) 

1 

1  +ae-*w 


1 _  =  1  =  1 

(1+ae  (l+ae^w)  1+ae  ^liJ+ae^aJ+a^  l+2aoos  (w) +ot^ 

The  final  result  is: 

Sy(w)  =  |H(eiw)  | 2  *= - -  (B.l) 

(1+a  )+2aoos(w) 

Autocorrelation  Function 


Starting  with  H(z)  for  this  case: 

1 


H(z)  = 


1  +az 


-1 


Per  Ref.  7:p.  158: 

=  h(n)  =  (-a)n*u(n)  for  a  <  l 

(u(n)  is  the  unit  step  function) 

Per  Ref.  5:pp.  391-395,  for  the  white  noise  input  case: 


Ry ( £ )  =  l  h  (n) *h (n-£) 

n=-°° 

Substituting  h(n)  above  into  Eg.  B.2 


(B.  2) 
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RV(M  = 


l  (-a)n* (-a) 
n=£ 


I  (-«) 


n-£ 


2n 


£  >  0 


£  >  0 


(-a)  n=£ 

The  summation  term  may  also  be  expressed  as: 

m 

,2n 


-  £-1 
2n  r  , 


l  (-a)  n  =  l  i-a.)  -  l  (-a) 

n=£  n=0  n=0 

Per  Ref .  12 :p.  8,  the  summation  terms  on  the  right  are  equal 
to: 


I  (-“) 

n=0 


2n 


- i— -y  (a  <1) 

1  -  (-a) 


£-1 

[  (-a) 

n=0 


2n 


1  -  (-a) 

1  -(-a)‘ 


2£ 


(a  <  1) 


As  a  result: 


l  ("<*) 

n=£ 


2n 


1  -  (-a) 


2£ 


(-a) 


2£ 


(a  <  1) 


1  -  (-a  )  1  -  (-a)  l-(-a) 

Substituting  and  using  (-a)2  =  a2  yields: 

2£  £ 

Rv(£  )  =  -A_  .  (~aL-  =  £  >  0  and  a  <  1 

Y  (-a)51,  1  -a  1  -a 


(B.3) 


1 


« 

5 

K 

Is 

a 


r. 

1 1 

j 

j 
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MAGNITUDE  MAGNITUDE  MAGNITUDE 


Power  Spectrum 


ALPHA=0.95 


ALPHA=0.5 


.16  0.00 
OMEGA 

ALPHA=— 0.95 


-3.16 


o.oo 

OMEGA 


ALPHA=0.8 


.16  0.00 
OMEGA 

ALPHA=0^5 


>3.16 


1.16  0.00 
OMEGA 

ALPHA=-0.25 


0.00 

OMEGA 


-3.16 


0.00 

OMEGA 


APPENDIX  D 


;***/,>  * 


ION  OF  THE  POWER  SPECTRUM  AND  AUTOCORRET.ATT 
FUNCTION  FOR  THE  FOUR  POLE  AUTOREGRESSIVE  MODEL 


Power  Spectrum 


|  H  (e3w)  |  2  =  H(ei“)  *H(e->) 


In  this  case: 


H(e_>)  = 


(1-ae  ^0e-^)  (l-ae^0e^w) 


,  i  (0-Ho)  -j(0-u)j  2  j2to 

1-ae  -ae  +a  e 


H(e3^)  = 


n  ~j0  “jw.  ,,  i0  -jto.  ,  -j  (0+to)  j  (0- to)  ,  2  -j2io 

(1-ae  e  ) (l-aeJ  e  J  )  1-ae  J  -aeJ  4a  e  J 


Multiplying  the  above  expressions  yields: 


1-oe3  (e-"'-ae-3  (fH“,«2e'32“-ae3  (0-H">+a2ej29+a2 

.  aV  (0-“).^-:  <6^>+a2«V326-aV3  (e-“> 


,  2  j2oo  3  j(0+co)  3  -j  (0+to)  ,  4 

+  a  eJ  -a  eJ  -a  e  J  +a 


Combining  terms: 


H(e3W)  ^(e"^  )  = 


l+2a2+a4-a3  (e'j  (0-w)+^  (6"w)  )+a2  (e^+e^20) 

-  a3(ej(0+w)^(0+w))4a2(e=>2u)+e^2w) 

-  a(ej  <0-u>)  )_a(ej  (9+to) ) 


Using  Euler's  relation  and  combining  terms: 


H(e3 w)  *H(e“jw )  = 


4  2  3  2 

14a  +2  (a  -cos  (0-0) )  [a  +a]+a  cos  (20) -cos  (0+to) 


3  2 

x  [a  +a]+a  cos  (2to) ) 


I 


Using  cos  (0- to)  +  cos(8+aj)  =  2cos  (  9  )  cos  (  w  )  ,  and  since 
Sy(w)  =  H(e^ w)  *  a2,  with  a2  =  1  the  final  result  is: 

sy(u)  =  — - 2 - 3 - - - Yl - 7 -  (D-1) 

1+a  +2  (a  -2  [a  -a]cos  (0)oos  (to)+a  (cos  (28)+oos  (2oj)  ) ) 

Autocorrelation  Function 


H(z)  = 


"  jlP1!  -j0  -1  2  -2  “  ,  j0  -10,  -1,2  -2 

l-aeJ  z  -ae  J  z  +a  z  1- (eJ  +e  J  ) z  +a  z 


Using  Euler's  relation: 


H(z)  =  - 1— _v  _2  |z|  >  |ot | 

l-2aoos(0)z  'L+a  z 

Per  Ref.  7:pp.  204-216,  partial  fraction  expansion  can  be 
used  to  find  the  inverse  Z  transform.  To  do  so  H(z)  can  be 
expressed  in  the  form: 


(1-ae  ^®z  (1-ae-^z  (z-ae  (z-ae^) 

Using  the  partial  fraction  expansion  and  table  look  up  [Ref. 
7:p.  158]  yields: 

7T 

h(n)  =  — : — t— r  oos(n0+0~4-)  -u(n)  for  a  <  1  (D.2) 

sin(8)  2 

Since  cos(e  -  1)  =  sin(e),  the  final  expression  for  h(n)  is: 

n 

h(n)  =  ,?  t  ,  sin( (n+1) 9) -u(n)  a  <  1  (D.3) 

sin (0) 
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k 

'/ 


For  simplicity  in  further  derivation  of  Ry(£)  based  on  h(n), 
Eg.  D. 2  will  be  used. 

Using  the  expression  for  the  autocorrelation  function  of 
a  random  process  represented  by  the  above  filter  with  a 
white  noise  input  [Ref.  4:pp.  391-395]: 


00  —  £  oo 

RyU)  =  l  h(n)  *h(n-£)  =  — - -  l  a  n-oos(n0+0-~) 

n=-<»  sin  (0)  n=£ 


cos( (n-£)  0+0~2) 


(D.4) 


n  =  £  in  the  summation  index  since  h(n)  is  causal.  £  is 
assumed  to  be  greater  than  zero  here.  For  £<  0,  Ry(£)  = 
Ry(-£)  by  symmetry  of  the  autocorrelation  function  [Ref. 
5:p.  388],  so  we  can  proceed  assuming  only  positive  values 


of  £. 


Using  the  trigonometric  identity  for  a  product  of 


cosines: 


-£  00 


Rv(£)  "  — ^ -  7  a2n[i  •  cos  (£0)+ioos  (2n0-£0+20-iT)  ]  £>0 

Y  sin2(0)  n=£  2  2 


-Z  °° 


— £  00 


R  (£)  =  °°sJ£e).,a  ..  .  £  G2n+__2_ -  l  a2n-  (2n0-£0+20-ir)  (°*5) 


y  2sin  (0)  n=£  2sin  (0)  n=£ 


£  >  0 

00  00  £- 1 

Using  \  ->  1  ~  1  and  standard  geometric  progression 

n=£  n=0  n=0 

identities  [Ref.  ]2:p.  8]: 


00  1 

v  n  1 

)  a  =  t - 

L  n  1-a 

n=0 


*--1  1  £ 

,  y  n  1-a 

and  l  a  3^“ 

n=0 


^Vvv'J  •  v  v  v  v  v  v  vjmr  ryyyjzs. 


For  the  first  term  in  Ry(£) 


2n 


n=!t 


00  J?  —  1 

r  2n  V  2n 
l  a  -  l  a 
n=0 


1-a 


21 


n=0 

21 

a 


.  2  :  2 

1-a  1-a 


1-a 


For  the  second  term  in  Ry ( l) : 


(D.6) 


let 


<J>  =  -  SL6  +  26  -  n 


Using  Euler's  relation: 


j  (2n0+4>)  -j  (2n0+4>) 
cos  (2n0+<J>)  =  - — * — - 


l  a2n.cos(2„9.«  -  t  I  a2V(2n8+'t"  + 

n=£  n=£ 


2  n =«, 


1  r  2n  -j  (2n0+<J>) 

+  J  l  a  ’e 

^  n=£ 


(D.7) 


For  large  n,  it  is  evident  that  the  a2n  term  will  tend 
to  make  the  term  in  each  sum  approach  0  for  a  <  1,  and  thus 
ensures  convergence  and  a  closed  form  expression  for  each 
term.  Pursuing  the  mathematics  required  to  find  this  closed 
form  expression  we  have: 


1 

2 


oo 


l 

n=l 


2n  j  (2n0+4>) 
a  *e 


J* 


?  (iej'')2n  =  l  (aej9)2n-Y(aej9)2n] 


l-l 


n=f. 


n=0 


n=0 
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1  r  2n 
T  1  a  * 


2n  j  (2n0+4>) 


j  4>  ,  ,  ,  j  6  N  2  £ 

ejy  r  1  l-(aeJ  ) 

2  jp  "  ~  JO 

1-  (oteJ  )  1-  (aeJ  ) 

a2£>ej  (20£+<|>) 

2  (1-  (ae-^9)2) 


(D.8) 


For  the  conjugate  term  we  must  have: 


,  00  _  .  2£  -j  (20S.+0) 

1  y  „2n  "J  (2n0+4>)  _  a  -e  JV 

2  L  a  *e  --in  2 

z  n=£  2(l-(cte  3V) 


(D.9) 


Next  a  common  denominator  must  be  found  to  sum  these  two 


terms: 


2£  j  (20£+<}>) 

U  »c 

2(l-(aej0)2) 


(l-(ae~j6)2) 

(l-(ae"j0)2) 


a2£ej  (20£+<)))_ct(2£+2)ej  (20£-20+$) 
2  (1-  (aej0) 2-  (ae~^8)  2-kx4) 


2£  -j(20£+4>)  „  ,  J6*  2. 

a  -e  J  ,  (l-(aeJ  )  ) 

2(l-(ae"j0)2)  (l-(aej0)2) 

2£  -j  (20£+4>)  (2£+2)  -j(20£-20+4>) 

a  e  _ -a _ e  _ 

2(l-(aej9)2-(ae':i0)2+a4) 

By  Euler's  relation: 


l-(ae^0)^- (ae  ^0)2+a4  =  l-2a2cos(20)+a^ 


Adding  the  terms  with  the  common  denominator  yields: 


a2£  [ej  (20£+4i)_a2ej  (20£-29+4>)+e-j  (2©£-»-0) _a2e~ j  (20£-20+<J>)  j 


noticing  that  the  same  result  must  hold  for  i  <  o  we  have: 


y*> 


^oos(lill9)  |oos((2-f-|il|)8-7T)-a2cos(|&|e-iT)  j  (D.12) 


2sin2(6)  1-a 


l-2a  oos(20)+a 
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GRAPHICAL  RESULTS  FOR  THE  POWER  SPECTRUM  AND  AUTOCORRELATION 
FUNCTION  FOR  THE  FOUR  POLE  AUTOREGRESSIVE  MODEL 
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Autocorrelation  Function 
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APPENDIX  F 


DERIVATION  OF  THE  POWER  SPECTRUM  AND  AUTOCORRELATION 
FUNCTION  FOR  THE  FOUR  POLE  AUTOREGRESSIVE  MODEL 
(WITH  TWO  POLES  ON  THE  REAL  AXIS! 


Power  Spectrum 

Using  initial  results  from  Appendix  D  with  the  necessary 
modifications  (including  8  =  0)  we  have: 


H(e-3W) 

H(e3w) 


Z  Taj  laT  i2u) 

1-a^e-1  -a^eJ  +otaa^eJ 


1-a^e  -^-a^e 


Therefore: 


H(e3u) *H(e”3u)  = 


~ia)  j-jw,  -jia)  J(jj,  2 
^ae  J  -a^d  J  -K^e  -aae 


2  -ja)  ja).  ,2  2  -ja) 

-  W  -%e  «aVV0aV 

,  j2a)  2  jai  2  ja).  2  2 

+  otaobeJ  -c^eJ 


Combining  terms  and  using  Euler's  relation: 


H(e3u)  *H(e“3  “)  = 


5  5  2 

1-  (cta+aj3‘Haaaj;)+aactb^  * 2oos  ^  +2aa0ibcos  ^2(jJ^  +aa 

+  ^aV^A 


Assuming  a2  =  1  and  since  Sy(co)  =  H(e3w)  •H(e“3aJ ) -a2,  the 
final  result  is: 


Sy(w)  = 


2  2  2 
1-2  (aa+at,+aaai:)+ota0ij;)) cos  ^  +2a_aj;oos  (2  u>)  +(^+20^0^ 


,  2^2  2 
+  Wb 


(F.l) 


Autocorrelation  Function 


H(z)  = 


~  j8-l  -j8  -1 .  -7 

l-aaeJ  z  -o^e  J  z  ta^a^z 


Letting  0=0 


H(z)  = 


1_  (aa+ab)  Z  1+aaV2 


z 

|z|  >  | a  I 


Expanding  in  terms  of  partial  fractions  we  have: 


H(z)  = 


z  +  % 


aa-ab  z'aa  Vaa  z_ab 


This  corresponds  to  the  impulse  response 


a  ol 

h(n)  s  ( — - —  *an  +  - *ar)  *u(n)  a  <  1 

aa”ab  a  Vaa  ^ 


(F.2) 


Proceeding  as  in  Appendix  D: 


Ry(£)  =  l  h(n)-h(n-£)  =  £  (— }  (_a - - } 

n=-oo  n=£  aa-ab  ab_a-'  a~-c^  a^_a- 
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■°b  V 
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03  a  a  _.ou~  an+l-£  n+1 
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Since  (a^)  =  -(o^-a  ) : 


(aa_ab)2  "  (Vaa)2  =  ~(Vab)  (Vaa)  =  “  (ab"aa)  (Vab) 


Ryd) 


- ~ — l  aln'a  i'1  l  <«  ^)n 

(aa_ab) 2  n=£  *  °  n=£  ^ 


-  Va"  J4(Va,n  +  ^  {£  1  0) 


Continuing  with  the  same  principles  and  assumptions  as  in 
Appendix  D,  we  have: 
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Making  the  appropriate  substitutions  in  the  expression  for 


Ry  (  &  )  / 

y*> 


we  have: 
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Combining  terms  yields  the  final  expression: 
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